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Abstract 


This thesis describes an investigation of the 


orientation dependence of the ae 


Na quadrupole relaxation in 
a Single crystal of NaNo, at room temperature. Nuclear 
quadrupole relaxation occurs via transitions where the 
Magnetic quantum number m changes by + 1 or + 2, the 
respective BY cpap iivttes of these allowed transitions being 
Wy and Wo- The ratio W./W, was measured as a function of 
the crystal orientation with respect to the external magnetic 
field using a steady-state nuclear magnetic resonance 
technique. 

The technique is an extension of the double-resonance 
method devised by Pound where the intensity changes of one 
of the quadrupole-split 23Na resonance components caused by 
severe simultaneous saturation of another component is 
measured. By measuring the intensity changes as a function 
of the spectrometer observing power, the effect of the 
observing power was corrected for using a linear extrapolation. 
In order to make accurate measurements for a range of crystal 
orientations, a sensitive flux-balance system was developed. 
This enabled the coupling between the spectrometer and the 
saturating signals to be eliminated. The final resonancés 


obtained with the equipment were such that their intensities 
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could be measured with an accuracy of 0.2% giving the. value 
of Wo/W, accune te, ato, 2. 0%. 

the wuadrupole relaxation, is characterized. by .a 
fourth rank tensor, the so-called M-tensor. Theoretical 
expressions for the orientation dependence of W/W, £o.1 
triad symmetry have been developed in terms of the M-tensor 
components by Pietila. Our experimental values of W/W, 
were fitted by a least squares technique to this equation 
resulting in a very good fit when the M-components in the 
equations for Wo and W, were considered equal. This implies 


Ui 


that the predominant relaxation mechanism is an indirect 


process. The results obtained were Mi 333/"1111 = 0.850 + 
Oe O55 M3333/M@1111 =-0.822 + 0.016, Mi473/1i11 = 0. L115 ae 
02028 “and My 3203/1411 =O .4477 410.017, wherer the: 2 taxis» is 


along the triad axis and the x axis is along the projection 
of a rhombohedral unit cell edge on a plane perpendicular 
to the triad axis. Using symmetry arguments it is shown 
that Mi 123 Should equal zero in NaNO... This is in marked 
contrast with our experimental value, indicating that at 
room temperature the crystal structures of NaNO. is mot as 
presently accepted. 


Previous measurements of the oy orientation 


dependence of => Na in monocrystalline NaNO, at 77°K by 
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Niemela using pulse techniques are reanalyzed. It is shown 
that agreement between experiment and theory is satisfactory 
and that the relaxation mechanism at 77°K is again via an 
indirect mechanism. 

In an attempt to interpret the experimental values 
of the M-tensor components, a calculation using a point 
charge model is presented. The calculated values of the 
ratios of the M-components are in surprisingly good agreement 


with experiment. 
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SECTION) 1 -iIntroduction 


Nuclear spin-lattice relaxation in solids occurs 
through the coupling of the magnetic dipole moment and the 
electric quadrupole moment with the lattice. For nuclei 
having a nuclear spin number I = 1/2, the quadrupole inter- 
actions are absent, but for I>1/2, both of these interactions 
contribute to the relaxation. In the magnetic case, the 
magnetic dipole is coupled to the fluctuating magnetic field 
at the nuclear site. In the electric quadrupole case, the 
electric quadrupole moment couples with the fluctuating electric 
field gradient at the nuclear site. 

The allowed nuclear spin transitions in the magnetic 
case are governed by the selection rule Am = +1 (Bloembergen, 
Purcell and Pound, 1948). In the quadrupole, case on the other 
hand, transitions corresponding to Am = +l and Am = +2 can 
occur Pound,” 1950) ; 

If the nucleus is situated in a crystal at a site 
possessing less than cubic symmetry, the time-average electric 
field gradient will be non-zero. This interacts with the 
electric quadrupole moment (for I>1/2) to modify the nuclear 
magnetic energy levels. In sodium salts the static quadrupole 


interactions of the ena nuclei, for which I = 3/2, are usually 
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much smaller than the Zeeman interaction of the nuclear magnetic 
moment with magnetic fields v104G generally available in the 
laboratory. This situation is usually referred to as the 'high- 
field' case and we shall be concerned with this case through- 
out this thesis. Since the static quadrupole interactions 
depend on the m value of the Zeeman levels, the 2I+1 levels 
become unequally spaced. For nuclei located in identical sites 
in a single crystal, the resonance spectrum will therefore 
consist of 2I components. 

Spin-lattice relaxation via. quadrupole interactions 
was first investigated by Pound (1950) who studied the 23Na 
spin system in a pure single crystal of NaNO, at room temper- 
ature. Pound showed theoretically that by severely saturating 
one component of the three=line resonance spectrum the in- 
tensity of other components would be significantly changed if 
the relaxation mechanism is quadrupolar. On the other hand, 
severly saturating one component would produce no effect on 
the intensity of the other components if the relaxation mech- 
anism is magnetic. By assuming that the quadrupole transition 
probabilities corresponding to Am = +1 and Am = +2 were equal, 
Pound showed that severe saturation of the central component 
of the spectrum (the so-called centre line) should cause an 


enhancement of the other components, the so-called satellites, 
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byjapt actor cot 43/2. +) Simileriy , ‘Pound ishewed that’ severe 
saturation of a satellite should depress the other satellite 
by a factor of 2/3 and enhance the centre line by a factor of 
Bife3's 

To check this, Pound performed a double resonance 
experiment in which one component of the spectrum was severely 
saturated and the intensities of the other components were 
measured using a marginal oscillator spectrometer (Pound and 
Knight, 1950, Watkins and Pound, 1951). Although the signal- 
to-noise ratios were poor, approximate agreement with the ex- 
pected intensity changes was obtained. The indication there- 
fore was that the dominant relaxation mechanism was quadrupolar. 
This was later confirmed by Andrew and Swanson (1957), by 
observing the saturation behaviour of the centre line and 
satellites. 

Abragam (1961) generalized the Pound theory taking into 
account the fact that the transition probabilities corresponding 
to Am = +1 and Am = +2 are not necessarily equal. Abragam 
pointed out that Pound's data implies that W,/W, is approximately 
1.5 rather than unity, where Wi and Wo are the respective 
probabilities of the transitions Am = +1 and +2 (Yosida and 
Moriya ji xl 95 Gide. 


In addition to assuming that Wa and Wo were equal, 
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Pound assumed that the effect of the spectrometer (observing) 
power was negligible. By taking into account the spectrometer 
power and lifting the restriction that W,=W., 


derived expressions for the double resonance behaviour of an 


Hughes (1966) 


I = 3/2 spin system. In particular, Hughes showed that the 
enhancement of a component should be linearly dependent on 
the spectrometer observing power when the observing power is 
small. It should therefore be possible to correct for the 
effect of the spectrometer observing power in the Pound 
experiment by measuring the enhancement as a function of the 
spectrometer observing power and using a linear extrapolation 
procedure. 

The first experimental measurement of W/W, in the 
high-field case was made by Goldburg (1959) for the 2393 spin 
system in monocrystalline NaNO, at 77°K. Whereas Pound 
measured the enhancement using a steady-state technique, Gold- 
burg employed a pulse method. The crystal was orientated such 
that its triad symmetry axis was directed along the externally 
applied magnetic field and W/W, was found™to “be 0N90"2770. 05. 

The aim of the work. discussed in this thesis is to 
investigate experimentally the orientation dependence of 
W/W, ina solid: “The investigation was Carried out on the 


23ua spin system in NaNO, at room “temperature. 
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A steady state method was chosen in preference to a 
pulse method (Hahn, 1949) because of the danger of interference 
with other components of the spectrum when the 90° pulse is 
applied to a given component. The quadrupole splitting is a 
function of the crystal orientation with respect to the exter- 
nal magnetic field and at certain orientations the resonance 
separation is so small that pulse methods can no longer be 
used. 

The double resonance method cannot be used for a gquad- 
rupole study of 23Na a9 NaNO, if the radio-frequency fields 
interfere with each other. We describe the design and con- 
struction of a sensitive flux-balance system that enabled the 
coupling between the radio-frequency fields to be reduced to 
very small values. If this had not been done, the measurements 
that we describe could not have been carried out. 

During the course of the work described in this thesis, 
the results of the orientation dependence of Wo/W, for Bares 
alee) NaNo. at 77°K were given by Niemela (1967) and compared to 
the theoretical dependence given by Pietila (1968). This work 
He ‘Critized in detail in Section iV. 

It was necessary for us to extend the theory in order 


to analyze our data and this is: done in Section Ir. . The 


equipment and experimental procedure is given in Section III. 
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Some auxiliary data concerned with the equipment is also given 
in Section III. The main body of the experimental data is 


presented and discussed in Section IV. 
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SECTION II Theory 


aiwn LRErOduction 


The application of a steady magnetic field Hy to a 
system of nuclei with nuclear spin number I and magnetic 


dipole moment wu gives rise to 21+1 nondegenerate Zeeman energy 


levels. These energy levels are given by 
En = ~MUH ) 
Hearts (Li 1) 
where m is the magnetic quantum number. In a state of ther- 


mal equilibrium, the ratio of the populations of adjacent 

levelcais, caven bystne, Boltzmann. factor exp | Ag | ieVenctse iV teas 
kT 

the energy separation of adjacent levels, k is the Boltzmann 


constant and T is the temperature of the spin system. For a 


Spin system at room temperature in an applied field of a few 
thousand gauss,-tiesboltzmann,factomis.of, the order of ort, 
So, to.an.excellent,.approximation..the. Boltzmann. factor can be 


written as l + cn ott follows theta Nuss the totals number 


of. spins, then. the.population of the mth energy, level is given 


by 


2 CLD. 2) 


eR Wile 1 og 2 OEE 
a ; bit F at ; 
De nag 
i 4 af ox 
j i oa 
a nota 

2 of. 44 bLesd ois ona vbeede 5 to ee 7, : 

Siionpem bas I eat alge weolouy Pere) isioum 2 

J 
| Yorsns Metieas stersiiopsiinoc [+1 lox ‘Selt bavie i some 9 
‘ yd. taevep exes ealavel eereAR and Loy 
_ a 
oe eis 4 < ~ 
|,’ 4S aR) | a aa 
n ar 0 State s ot .2edmur MUS TEND | i isbeeuilea ater scale 2 
an A 
, Pa’, dnsostbs to ancitsiyqog ons to, oigat ous: » ies. : 


> | seh fA siodw te | axe xotdst nntems#kOS erid uch aowie, a | 
4 Ea _ 
ie) nein L6G sft Bi gi palevol inode ee. to otgnasase xt " rus 3 


1g 8 40% meteye nigqe Sds to citi si py od pe ‘ba. 


ibs | Lpet te ses ae 


NMR occurs when photons of frequency v given by the expression 
E/h interact with the spin system. A quantum of energy will 
exite transitions between the energy levels if it has the 

Same valuetasPthepleved spacings. 

Let us consider an alternative picture. The appli- 
cation of a steady magnetic field will result in a precess- 
ional motion of the spin about the direction of the applied 
magnetic field. Such a motion will be characterized by the 


Larmor frequency of precession that is given by the equation 
Qi) =SythH CLIZS)) 


where se is the nuclear gyromagnetic ratio for the particular 
Species. Let-us suppose; thatia circularly polarized field is 
applied such that its plane is perpendicular to the steady 
magnetic field Ho: The frequency of the polarized field can 
be adjusted until the rotation rate of the radial vector 
representing the field as equal) tothe precessional rate of 
the spin about the field Hy and in the correct sense. In 
this case there will be a sustained interaction between the 
spin. and the circularly polarized field resulting in the re- 
orientation of the spin into another energy-state. such an 
interaction is the basis of the Nuclear Magnetic Resonance 


phenomenon and this description is equivalent to the energy 
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level picture given in the preceding paragraph. 

The Magnetic transitions resulting from absorption 
of energy by the spin system and resulting in the. redistribu- 
tion of the spin population is governed by the selection 
Piles’ Meo nt ol -eThe sprobabilityeforssuchaa trensdtion,takes 
the form (Bloembergen, Purcell and Pound, 1948 henceforth 


referred to as BPP, 1948) 


2 


P =, Pp = (3) ‘é 


Hf g(v) [I+m] [I-m+1] Cee -Ay) 
m?>m=— 1 m-1l>m 


where g(v) is the normalized lineshape function and Hy is the 
Magnitude. of the circularly-polarized applied field. The net 
absorption of energy by the spin system occurs because of the 
slight excess of population in the lower energy levels and 
the intensity of the nuclear absorption is directly propor- 
tional to the. population difference. between the two levels 
between which the spin transition takes place. 

For an isolated nuclear spin system, devoid of inter- 
actions with another system, the continued presence of the 
appropriate radiation would in time lead to equal population 
of the Zeeman levels and there would be no further net energy 
absorption by» the spin system.. In such,a state, the charac- 


teristic spin temperature will be infinite and the spin system 


is described as being completely saturated. 
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In practice the spin system is not isolated but in 
fact interacts with the lattice. Such an interaction, the 
so-called=spinslattice interaction, has a finite probability 
of inducing transitions. However, in this case the probabil- 
ity of downward transitions exceeds that of upward transitions 
(Andrew, 1955). " yClearly. such “spin-lattice interactions will 
form the basis of the relaxation mechanisms for the particular 
nuclear spin system. 

To illustrate this let us consider the simple case of 
a spin system with I = 1/2. Using equation (II 4) we write 


the transition probability as 
- ~ 
4 n (ele 5) 


Defining the upward transition probability arising from the 
spin-lattice interaction as W we can write the downward tran- 
Sition probability as W(1+A) (Andrew, 1955) where A = hwkT, 
hv being the energy separation between the two levels. By 
probability arguments we can then write rate equations for 


the populations of the two levels as 


No-= N, W(1itA) + NjP, = Ny Wit= Ny xe 


i ui i i i 
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wnere the N’s are, the populations of the particular Pevyels. 


Since sA<<1 and Ny = Ny = NASSN; 1t) follows’ that 
ee 
= ad + a 
Ny WN , Ny ene (tr 8) 
2 
and 
pil a eae oa epics Oe. (EE. o) 
Z 
Hence, it follows that 
N, = ~2WN, - 2Ny Woe eoNee 6 (II 10) 


Now uSing equation (II 2), NjA is the equilibrium value of 
N, which we denote as n and therefore 


Be pest (nQ7 Ny) _ one Cr i es 


For P =0 the Solution for this type of equation is of the form 
fe) 


Nj=an,[i-c exp (-2Wt)] where C is a constant that is determined 
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by the initial conditions. The time constant for the ex- 


pPonential return ytovequilibrium is Crh which is defined as 


the spin-lattice relazation time T,- 
inMmMeneypnesence: Ol arradieation field, att. forrows rom 


equation (II 4) that the steady-state solution of equation 
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If the applied radiation field is sufficiently large such 


that Ye Hf 


severely saturated. This derivation of equation (II 12) was 


TS g(v)>>1l then the spin system is said to be 


first given by BPP (1948). However, for NMR, in solids 
equation, (21 22) 1s only valid: for small values of Hy except 
at the centre of the resonance where it is true for large 
Hy- This was demonstrated by Goldburg (1961) for ona in 
NaCl. 

For the case I = 1/2, spin-lattice relaxation is due 
to the magnetic dipole moment of the nucleus interacting with 
the time-dependent magnetic field produced at the nuclear site. 
However, measurements by BPP (1948) indicated another type of 


relaxation process. In water doped with deuterium, T) for 


a proton is approximately 3 seconds whereas Ti for a deuteron 
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is approximately 0.5 seconds. Since the magnetic dipole for 
Ppvesieuceron 1c much. less than thet of the,proton,, 16 ais 
elear that another relaxation mechanism,is involved. . The 
deuteron has an electric quadrupole moment and it was pro- 
posed by BPP (1948) that it was this interacting with the 
time-dependent electric field gradient that caused the short 


ak In the case of heavier nuclei which have large quad- 


Ll: 
rupole moments, broad linewidths have been found that indi- 
cate relaxation times as short as 107 *seconds In uguLds:, 

In the case of liquids, the quadrupole interaction, 
although broadening the resonance line, does not result in 
any additional {structure ‘of the: line. Such» structural changes 
would arise from the interaction of the quadrupole moment with 
the time-average electric field gradient in the crystal. How- 
ever, in the liquid the local configuration time is suffi- 
cently short that such an interaction-is averaged out. In 
the case of a solid, however, the local configuration 
does not fluctuate so rapidly, and a non-zero time-average 
electricitiretd gradient may existvat the nucilearesi1ts.. Thus, 
the interaction of the electric quadrupole moment with this 
field gradient may result in fine structure of _the resonance 
(Pound, 1950). In passing, we mention that the time-averaged 


electric field gradient at the nuclear site is zero if the 
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HUCKFeUS™IS in a Site possessing cubic symmetry. On the other 
hand electric quadrupole relaxation may still occur. 

Let us now consider in detail the quadrupole inter- 
action in a solid and the subsequent effect that such an inter- 
action will have on the Zeeman energy levels arising from the 
application of the steady field Ho to the system. Consider 
a nucleus with atomic number Z whose total electric charge is 
distributed over the nuclear volume with a density 0 (r). 

Let V(r) be the electric potential arising from all the other 
charges in the environment excluding that of the nuclei 
under consideration. The electrostatic interaction H is given, 


over the nuclear volume, by 
H = fo(r) V (r) dv. Garni) 


The term V(r) can be expanded over the nuclear volume about 
the nuclear centre of mass by a Taylor series giving 


:: oer er 
H = fd-xo(x) |"o* Zs Ee ae + 1 bay | 2 pects | (pti 4) 
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>a a 
where the components of x in the Cartesian representation are 


1 2=Y and X32 with the summation over each of the sub- 


Scripcs,extending, over 1,2 and 3. Thus one. obtains the ¢x- 


Xo=XGe x 
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pression 
" 2 
H = ZeV +. ¥4EPY [ww] +22, Q. pee Bo aie Cia sr) 
2 J J lgxto 3 aaa Oh. OO 
J Speier 
where 
she 7 
oe = fd BD gst) Xs is the nuclear electric dipole moment and 


(£2516) 


Qo = fa°xp (x) x %, is the nuclear electric quadrupole moment. 
(araheealay 7s) 


The first term in equation (II 15) represents the electrostatic 
energy of the point charge nucleus and since it is independent 
of the nuclear shape or orientation it is of no interest in 
this work. The electric dipole interaction term vanishes by 
virtue of a parity argument. The basis of the argument is 


that the ground state wave function ee has a definite parity 


such that ¥,. Panay is unchanged if all nuclear coordinate signs 
; a — 
are reversed. Thus, it follows that RN =? 9 ance ay Cry nit 


Since the term ea in equation (II 16) has different signs in 
the different octants then the term as vanishes. This 
argument also holds true for the case of the electric octupole 
term so that in fact the next non-zero term in the expansion 


is the electric hexadecapole term. One can therefore write 
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the interaction Hamiltonian in the form 


A 


Ho : Be Q-% Mite + hexadecapole terms +---- (ii 13) 
where 
2 
ie rae 
jk Sp ae wc ] 
Hage 


Looking at the approximate orders of magnitude of the 
interaction terms in equation (II 15), the electrostatic energy 
Ze*/x, is MAD PonTS having taken ry to be typically one Angstrom. 
The quadrupole term is ver“ (e/r,) = eV (x /r.)* where roe 
the nuclear radius is v10+ems. It follows that the quadru- 
pole interaction B10) of the monopole interaction. The 
hexadecapole interaction is some ome of the quadrupole inter- 
ackEion "ruling it, out for all intents and purposes. 

We shall henceforth be concerned with the so-called 
"high field" case where the quadrupole interactions are small 
compared to the Zeeman interaction. The effect of the quad- 
rupole interactions can therefore be treated by perturbation 
theory. A first-order perturbation treatment was carried out 


by, Pound: (1950) and ‘results in the nuclear energy levels being 


given by the equation 
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— 2 ‘ == 
‘it eG be = (EE) | Voz. (Tie) 


I 4T (2I-1) 


where eQ is the nuclear electric quadrupole moment and ve 


is the time-average zz component of the field gradient at the 
nuclear site and zis the direction of the applied magnetic 
field. In the absence of quadrupole interactions, the Zeeman 
splitting produces 2I+l equally-spaced energy levels. As can 
be seen from equation (II 19), quadrupole perturbation results 
in a degree of perturbation that depends on the magnetic 
quantum number m of the nuclear level with the result that 
the energy levels are no longer equally spaced. If the envi- 
ronment has axial symmetry then the field gradient term is 


given by (Pound 1950) 


Vere = [2] eg [3 cos* 8-1] Cr) 0) 


where 6 is the angle between the symmetry axis of the crystal 
and the external field Ho: Also, eg is.a scalar descriptive 


of the electric environment defined as 


eq = fal3cos~ oF | au dav CEieces) 


where the integral is taken over all the charges outside of 


the nucleus, ris the vector joining the nucleus to the 
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volume element dV, o is the charge density and ¢ is the 
angle between r and the symmetry axis. From equation (II 19) 


it follows that the component lines have frequencies 


Vv = + 3e"0g (2m-1) (cee oon) 


Pencil Brersinn (II 22) 


where ve is the unperturbed Zeeman frequency. 


For a single crystal of NaNO,, all the 22a nuclei 
are in environments possessing axial symmetry and identical 
23 


Van. Since. [aay 2’ for the 


uh Na nucleus, and since the symme- 


try of the crystal is lower than cubic, the tena nuclear 
magnetic resonance will consist of three lines. These are the 


so-Calleéd centre lihe arising from the” transition m=1>"-1”"and 


Piast Ly 
the satellites corresponding to the transitions m=3 > 1 and 
2a tee 
Mgeo-S ~ -l. he tirst-Order perturbation treatment! of “Pound 
27 2 


has been extended to higher orders by Bersohn (1952). The 
frequency v of the centre. line and the frequencies v, of the 
Cc + 


satellites are given by 
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where p is the dimensionless parameter e“Qq/hy,. The first- 
order term in 9 give the equations as expressed by Pound 
(1950). To this approximation, quadrupole interactions affect 
only the satellite frequencies, with each satellite being 
affected the same degree. Second-order and higher even-order 
perturbation terms affect the centre line and also the 
Satellites shifting the satellites with respect to the fre- 
quency Vv. whilst maintaining a constant satellite separation. 
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For Na in NaNO the third-order. term is: of the. order, of 


td 
10Hz and for most purposes can be ignored. 

The frequency separation of the satellite components 
as a function of the orientation of the symmetry axis of the 


crystal with respect to the direction of the applied field 


Hoe is to first order given by 


vWgtarNtiog helteg7eq 3cos* 6 -1 Poser 
Z hy, | 


where the quantity e“Qgq/h is usually referred to as the 


quadrupole coupling constant. 


b) Population dynamics for a, spin system with I=3/2 


We consider a system of nuclear spins with I=3/2 
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Situated in a single crystal on sites which have identical 
and non-zero Vis. . The resonance spectrum therefore consists 
of three lines. Furthermore, we shall suppose that the re- 
laxation of the spin system is due to homogeneous quadrupole 
interaction. The allowed spin transitions are governed by 
the selection rules Am=+l and Am=+2. 

We shall suppose further that the spin system is 
subject to a double resonance experiment of the type carried 
out by Pound (1950) in which one resonance is severly satu- 
rated whilst another is simultaneously observed using a much 
weaker Signal. In his experiments Pound investigated the 
aoNa spin system in NaNO, at room temperature and showed that 
the relaxation mechanism was predominantly quadrupolar. This 
was later confirmed by Andrew and Swanson (1960) who studied 
the saturation behaviour of various components of the spec- 
trum. By assuming the probability of the two possible quad- 
rupole transitions were equal, and that the observing power 
of the spectrometer measuring the component undergoing the 
intensity variation was. having no. effect on. the spin system, 
Pound showed that the intensity of the satellite should be 
enhanced by ai factor. of 3/2 when the centre line transition is 
severelysaturated. Similarly, Pound showed that when a sat- 


ellite is severly: saturated the.centre line should be en= 
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hanced by 5/3 and the other satellite depressed by a factor 
of 3/2. Experimentally, Pound obtained enhancement values in 
rough agreement with these values. 

We now generalize the Pound theory to take into 
account other considerations which are necessary to carry out 
accurate double resonance experiments. The quadrupole in- 
duced spin-lattice transitions probabilities are given by 


(Yosida and Moriya, 1956) 
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ae : (reoG) 
2I (2I-1) 


The downward quadrupole transition probability by the same 
selection rules are given by the above equations multiplied 
by the appropriate Boltzman factor. 

From examination of- the unit cell of NaNO ,1it is 
apparent that- for relaxation purposes there are in general 
two different Na sites though they both have the same V__ 
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and therefore the same quadrupole splitting. The quadrupole 
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Figure II a shows the nuclear transitions considered 
in the setting up of the population rate equations for the 
I = 3/2 spin system. Transitions resulting from the inter- 
action of the electric quadrupole moment with the time- 
dependent electric field gradient are labelled in terms of 
the quadrupole transition probabilities Wi and Wo transi- 
tions resulting from the interaction of the magnetic dipole 


moment with the time dependent magnetic field are labelled 


in terms magnetic transition probability W3- 


m <sinks.. ae ? 
Se. seers _ pelladst S18) Pees nieie aid 
“Lens 13 ow Bris ” ME emo 
a | he ama s16 bisLt aitomen eo vip 


ow wit ada noitianens + cits 


W[I+A(1+8) ] 


23 


transition probabilities Wy and We given by equations (II 25) 
and (II 26) are therefore in general different for both sites 
and our theory must take this into account. 

We shall consider the general case where Wy and Wo are 
moequal. We shall further introduce into our, calculations, Gt 
the population dynamics the refinement that the fia energy 
levels are not equally spaced. Ignoring second-order frequency 
shifts, we express the frequencies of the components of the 
ena spectrum in monocrystalline NaNO. as Vv) for the centre 
line, Vv, (1+) for the high frequency satellite and v, (1-9) 
for the low frequency satellite. The two spin systems, say 


A and B, will be coupled by spin-exchange interactions since 


; 23 
corresponsing Na resonance frequencies are equal (Andrew 


and Swanson, 1960). We label the intrinsic spin-exchange 
transition probabilities as @ for m = -3°-1,8 for the tran- 
aaa 
Sition a = -1> l.and.y form m.=_3>,1. .We now develop the rate 
Ze 2 ee 


equations for the various levels of the A spin system, and 
calculate the enhancement of the high frequency satellite 
where a saturating power > is applied to the centre line, > 
being given by Lye Hig(v). This will correspond. to the ex- 


perimental case where the observing power of the spectrometer 


is zero. The allowed transitions are shown schematically in 


Figure II a. 
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i) m= -3/2 


Using Nae for the population of the m = -3/2 level 
of the A spin system and similar labelling for the other 


levels, the time-dependent population of the m = -3/2 level is 


Se eA, A Lg z 
WZ = Ny Wy + Np WS Nz Wy) (1 + A) Nz Wo (1 + 2A) 
2 2 2 2 


This can be rewritten (c.f. Andrew and Swanson. 1960) 
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sik 
Ny = Ng Wy (1 + A) +.N3 Wy - Np WL - Ny WS (1+ 2A) 
2 3 2 2 2 


sg 0 a a Nr - NG + 8 Ne - Ne 
4 4 ° (TE 28) 


Similarly Ny is given by 
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Wy (1 + 2A) - Ny wh (1 + A) 
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(Eise30) 


Using, equauions (FT 427,),4 (BL 28)*, +(01 429) mands (Il 30))wexget 


following expressions for the population differences 


NU = NO (4P*) + ND (-2wy-wWo) + NS (WO) + an’ (2W4) 
aoe 
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A , A 
Ne = NO (4P*) + NY (WO) + ND (-2wi-wO) + an® (aw) 
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We wish to find the steady state solution where Ny = Nia 
we = 0. The term P~ can be eliminated from these expressions 
£0 Ssimplary obtaining the solution= = Multiplying Aili) by 2 


and adding “to (11,33)) wevget 
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bakewise multiplying Cll 32) tby 2 amd-adding to, (11 33) we 


have the form 


= (ey. 
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) +N (3a) + 2n, [wit 


é (w+). (hr 35) 


A similar set of rate equations will exist for the spin 


system B and by similar manipulation they can be put into 


the form 
0 = ND (-2wS) + Ne (-3Wi-3W5-3y) + ND (WS+We +a) 
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If we assume that the power P~ applied to the centre line is 


sufficiently large to cause complete saturation, then we can 


set the term Ne = Ne = 0 thereby giving four steady-state 
equations: 
Gene eu Te chy ee coluae ame ToPTa te) octamer) 
a ii 2 A = ale: Hh 7 1 a 
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We now solve these equations in order to obtain the 


value no and n® which will be directly proportional to the 
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Spin systems A and B. The total intensity of the enhanced 


at A B 
1? Ehe_ sum of Ny and Ny - 


Assuming a,8,y>> We Wee Wor Ws, the solution takes the form 


satellite component is written as Ne 
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Since the unsaturated intensity is given by 
Nr = 2n. (ite), acim eHep) 
the enhancement is given by 
Es [cw ten) + (We+2WS) + S(WhHWS) + $ (WW) 
(1+6) [wy + We + We + Wo] (ii 741) 


It is now convenient to introduce mean transition probabilities 


Wi and Wo defined by 
(peas he and WwW, = We + we 
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When 6 = 0, this reduces to the single spin expression 
(Abragam 1961, Hughes, 1966). We can therefore treat the 
23 


Na spin system in NaNO, as a Single spin species with tran- 
sition probabilities Wi and Wo- 

Let us. now.develop the equations .for a single spin 
species allowing for unequal energy level spacings and as- 
suming that radio-frequency power is applied to all three 
resonances simultaneously (Hughes, 1966). In this case the 
Spin=exchange interactions need not be considered since we 
exclude the crystal orientation where the resonances overlap. 


The rate equations for the various levels now take 


the form: 
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where Pay Pi and Pay is the power applied to the Nj, N, and 


N_, resonances respectively. It then follows that 


Ny =lte (4P °) & Ni (-6P) - 2W,-W.) + N_y (W.) 
+ AN (2W,) + ANG (2W, + 2W.), (II 48) 
om i —— a: wy) 
4 4 
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We now normalize the rate equations with respect to 


Wie putting the ratio Wo/W, = X, Pi /Wy = oe, Pi/W, =P, and 


P’,/W, = P_,- Thus the equations take the form, 
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N_)/W, mi Neeley NW) + NL) (=6P.,-2W Ww.) +n, (2W,) 
- no (—2W, a 2W.) Cry 52) 

and 

N/M = N, (-8P,- 2W,) + N, (3P,+W,-W,) + N_, (3P_,+W,-W,) 
tpadh (4W., = 2W,) « (Tele 53)) 


Let us now consider the situation where P =P_j=0. 
Such a situation covers the experimental case where the high 


frequency satellite is observed using the spectrometer. The 


population difference Nera en which determines the intensity 


of the N, resonance is given by 


dh 1 


(II 54) 
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Let us now consider the situation where, in addition, the 
centre line transition has power applied to it such that 


Po?>l, P This will correspond to the double resonance 
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situation where the high frequency satellite is being ob- 
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enhanced intensity of the N, resonance is given by 


= 


1) (rr 55) 


We Po] (2x + 1) + 6(1+x) ] (1+x+3P 


the enhancement BE of; the high frequency satellite is therefore 


given by 


E = (2x+1) + 6 (1+x) » | [ 22430 + Po 3 (2+x) ] 


| (II 56) 
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Expanding equation. (11 56) to first order ain Py as allowed 


£Oxr Pi<<l corresponding to low observing powers of the 

spectrometer, the enhancement becomes 

BE =< }1+2x 1+3(2+x)P) - 6x 1 + 3P, zl CEEVS7) 
1+x 2 (1+x) 1+2x Ere 


Similarly, the enhancement of the low frequency 
Satellite as a result of severe saturation of the centre Line 


transition can obviously be expressed to first order in the 


power Pia as 


oul 
E = ees 1 + 3(2+x) P_y Lit 6x 1+ 3P_, 
1+x 2 (1+x) 1+2x =e (Gast 58) 


Similar expressions can be obtained for the enhance- 
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ment of the various components due to the saturation of any 
other component. However, those of interest to us are the 
two given above in equations (II 57) and (Il 58) “and a ‘third 
case corresponding to the enhancement of the centreline due 
to severe saturation of the. high.frequency satellite. The 
enhancement for this case expanded to first order in the 


observing power is given by 
Ber ths+2xebeldabe2ileeBe |/1 + S(iex)||2 4+ atin) |". 6a, 
2+X x 3+2x x (2+X) | 6 


For simplicity we shall henceforth refer to the three en- 


hancements.givenpby.(Lin5/),, (11.58), ,ande(Gip59)--sassGase I, 
Case II, and Case III respectively. 

In all three cases the enhancement, for low values 
of the observing power, has a linear dependence on the 
observing power. In the limit of zero observing power, the 


enhancements for the three cases,I,II, and III are respectively 


Ev= 1+2x 1+ 6x CET 60) 
Lax Lt 2285 
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1+x 1+2x 
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Thevalue 6 can be found from measuring the resonance frequen- 
cies using the relationship ii la 2v 05. If the enhance- 
ment at zero power is found by extrapolating the linear portion 
of the enhancement/observing power curve, the value of W,/W, = Oe 
Cone Desroundslromscquatons, (Il 60), > (LE 61) and. (LL.62):. 

This is the basis of our proposed double resonance 
experiment, namely to measure the enhancement of a component 
due to severe saturation of another component as a function 
of the spectrometer power. 

It has been pointed out by Snyder and Hughes (1970) 
ehatsfor the case: ot quadrupole interaction of yan indirect 
nature, such as harmonic and anharmonic Raman processes, 

We ey and Neat are independent of m. If, however, the re- 
laxation is caused by a direct process such a one-phonon in- 
teraction then Ne ea and Wier eo will depend on m. We, 
therefore consider the possibility that Wy and indeed Wo will 


be different for the different transitions. Using a notation 


, : 2 
in keeping with that developed previously, we define Wy LOT 


the transition m = -3+-l and W. for the case m = -3+1 with W 
a 2 = 5 al 
" Zee 2=2 
and W., for the transitions m = 0-1 and mm =,3>-l respectively. 
22 2502 


Further we define the relationship between them as 


W, =, (1+) and w= W, (1+8) (It 63) 
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where 8 and > are small quantities. The appropriate rate 
equations can be set up in the manner described earlier. To 
first order in Poy the enhancement of the high frequency 
Satellite caused by severe saturation of the centre line 


transition has~the ‘form 


2 i = on 7 
CS 2W.+W, 1 W,W.¢ (SW, 2WoW,) 8 | W569 
W,+W. (Wj +2W.) (W,+W.) 2(W,+W.) (W,+2W,) (2W.+W, ) 
2 

MLth.+ P,3(W.+2W,) 1 - (Wi +2W,) > ~ W,W.8 

2W. (W.+W. ) 2 (W.+W.) eis 

Pvp cag hee 2(W,+W,) (Wo+2WW, 
-l 
Bo12 Le [ee 


For: the direct phonon interaction mentioned above, the phonon 
density of states should increase monotonically with frequency 
and therefore W the quadrupole transition probability for the 
high frequency satellite should be greater than W, the quad- 
rupole transition: probability tor the low frequency ssare lice. 
The quantities 6 and.¢ defined-in equation (II 63) are there=- 
fore positive. It can. be seen from (ITI 64) that the etfect 

of these quantities for the situation Pi= 0 is to reduce the 


enhancement of the high frequency satellite. By symmetry 
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arguments the enhancement of the low frequency satellite 
due to severe saturation of the centre line can be obtained 
£rom equation (il, 64) by changing the’ signs of 6, @ and 6: 
The effect of the quantities 9 and ¢ in this case increases 
the extrapolated enhancement of the low frequency satellite. 
Finally, we consider the correction that should be 
applied to the measured enhancement values to take account 
of the fact that the saturating power is not in practice in- 
finite. (We shall see later that such a correction is 
necessary to obtain the required experimental accuracy). For 
the case were the centre line is severely saturated witha 
power Py it can be shown from the expressions presented by 
Hughes (1966) that the saturated enhancement of a satellite 


is given by 


ae es (II 65) 
O x 


in the limit of zero observing power and ford=0. By 
measuring E as a function of the saturating observing power 
Ai, with a very small power, and fitting the data to 

an equation of the form given by (II 65) a correction for the 


non-infinite saturating power can be applied. So long as Po?rl 
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the error in such a correction factor due to the non-zero 


observing power used in the measurements is small. 


c) Quadrupole spin-lattice interaction 


The Hamiltonian for the quadrupole spin-lattice in- 


teraction can be written as 


L 7 Q (II 66) 


where He is concerned with the lattice states and Hy is. the 


zeeman term. The Hy term represents the interaction of the 
nuclear quadrupole moment with the electric field gradient 
produced by the lattice. The electric field gradient component 


V at the nuclear site may be written as 


ik 


= as ’ I 
cone a + et (oes Os7) 


The term ae is the time-average field gradient whose inter- 
action with the quadrupole moment determines the static quad- 
rupole interaction. This interaction, as we discussed earlier, 
Wil. result in,a,perturbation .of the, Zeeman, energy levels 


resulting in unequally spaced energy levels giving the compo- 


nent structure of the spectrum. The term Mia, WE) is the time- 
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dependent field gradient whose interaction with the quadrupole 
moment will determine the quadrupole relaxation. 
Following the notation of Pietila (1968),we express 


the quadrupole interaction Hamiltonian in the form 


oe ane acd ba i 
2 y=-2 (II 68) 
where 
By a, le 
2 
ln 
B = th 6:) wae # ee 
ee Cli/v.6)) dak 367 = Viel, 2511 VAseu 
a erent a 4 evan 
fe) 2 
On = k[3I0 Se Oe ee 
(II 69) 
ob st (y.67 20k CL ay + rE he 
Q°* = (¥6/2)kI%, 
and k = eQ/[21I(21I-1)]. 


The term gr involves the nuclear quadrupole moment and the 
nuclear spin parameters I, and I,. The subscript yu takes the 


five values 2, 1, 0, -1, -2 representing the five possible 


y i) * ve - x 
“ aeeageo ewig(BOCI) skisota col pre a ania els 
| mxot osit Mi asin? 2 tme sostoszedat : ABs 


(83 Tt) Ss ol | 


: Pax” 
‘sae aol & 


; Hy” » tea hs ale my i 


(ben) - Sena = i: 


= eee ae 


(@9 11) | | 
‘ Gry tlds epee » oe 

. TA(SNO\) = ig a 4 
| j 


> 4 | r al 
{firs von a * | 


iS | San 


4l 


Caption Ji 


Figure II b shows the Feynman diagrams of the relax- 
ation processes. The so-called direct process is shown in 
i), the indirect harmonic Raman process in ii), and the 
indirect anharmonic Raman process is shown in iv). The 
process indicated in iii) is a three-phonon process which, 
as is discussed in the text, is not an effective relaxation 


process. 
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Soin eranstertons vam t2, A m=) el Gand A ms=.0.- The: term 
Br PepLesences thecelectricetveld. gradientrand this can be 
expanded as a power series in terms of the displacement, 
relative to the relaxing nucleus, of the ions representing 


the lattice. Such a power series has the form 


AB Pe ee ee eee (Tr 70) 


where La represents a component of the relative displacement 
from the equilibrium configuration. 

It is convenient to categorize relaxation processes 
according to whether they involve the displacement Ch linearly 
via the second term in equation (II 70) or quadratically 
vra the third ’term in (IL 70); etc. (Van Kranendonk, 1954). 

Ns) an “examplie of the. first category is the directiphonon 
process in which each spin transition is associated with the 
creation or annihilation of a lattice phonon (Van Kranendonk, 
1954). Such an interaction is represented by a Feynman 
diagram in Figure II bi. However, this process should be 
relatively ineffective (except at low temperatures), since 
the density of phonon states with frequencies equal to the 


nuclear resonance frequency (v10/Hz) is small. 
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As an example of the second category is the indirect 
harmonic Raman process first discussed by Van Kranendonk (1954). 
The treatment consisted of a first-order perturbation calcu- 
lation involving the interaction of the incoming phonon with 
a nuclear spin to produce a spin transition and an outgoing 
phonon of different energy. The energy difference between the 
incoming and outgoing phonons is equal to the energy difference 
between the initial and final spin states. This process is 
represented by the Feyman diagram in Figure II b ii. In this 
process all phonons are involved thereby making the process 
in general much more effective than the direct phonon process 
BeLterredr LO, above. 

Van Kranendonk and Walker (1967,1968) have proposed 
another relaxation process that they have called the anharmonic 
Raman process. The total potential energy of a crystal can be 
expanded aS a Taylor series in powers of the atomic dis- 
placements from the equilibrium configuration. The quadratic 
term in such a series is the harmonic description of the 
lattice and higher powers are called anharmonic terms. The 
first ankarmontc (eubile)term “gives rise (to the three—-phonon 


process, shown in Figure II b iii, which plays a role in 
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lattice thermal conductivity. The anharmonic process con- 
Sidered by Van Kranendonk and Walker (1967, 1968) is a second- 
order process involving the first anharmonic term in combi- 
nation with the linear term in equation (II 70). Sucha 
process therefore falls in our first category. The Feynman 
diagram for such a process is shown in Figure II b iv. The 
phonon causing the spin transition in the anharmonic Raman 
process is a virtual phonon as opposed to the real phonons 
involved in the lattice conductivity process shown in Figure 
Pio. thes third: real pnonon: of Figure «11 bili was pro— 
posed as a Suitable phonon for spin relaxation by Kutsishvili 
(1967). However, Van Kranendonk and Walker (1967) have pointed 
out that such a process is quite ineffective as a relaxation 
mechanism. 

The temperature dependence of the spin-lattice relax- 
ation time T) associated with indirect harmonic Raman proc- 
esses and anharmonic Raman processes is the same as for the 
magnetic relaxation processes due to lattice vibrations. If 


one assumes a Debye frequency spectrum, then the temperature 


dependence is 
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ATs for Peed 0102 5 (it a7) 


where a and b are constants and 85 is the Debye temperature. 
Measurements of Ty (Wikner, Blumberg and Hahn, 1960, Weber 
and Allen, 1962) have confirmed the r dependence of Tie 


However, the p/ dependence has not been verified because 


of magnetic relaxation process at the low temperatures. 


d) Orientation dependence of the quadrupole spin-lattice 


relaxation 


As is seen from equation (II 69), the BY and hence 


the quadrupole relaxation probabilities depend on ey Vee, 


XX 


ee etc., where the z axis is along the magnetic field di- 
rection. (A coordinate system where the z axis is set along 
the direction of the externally applied field is later referred 
to as the K~ coordinate system). If the crystal is rotated 
with respect to the magnetic field, these field gradient com- 


ponents will be different and hence the relaxation probabil- 


ities Wi and W. should be orientation dependent. 


The orientation dependence of the transition prob- 


abilities was first considered by Pietila (1968), for a harmonic 
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lattice. Pietila considered the quadrupole relaxation due 
to the direct process and the indirect harmonic Raman process, 
for which he derived the orientation dependence for 3-fold, 
4-fold and 6-fold crystal symmetry. The question arises as 
to whether the form of the orientation dependence for the 
anharmonic Raman process is the same as that calculated by 
Pietila. This question was considered by Snyder and Hughes 
(1970), who showed that the orientation form was independent 
of the relaxation process. The following discussion follows 
the notation used by Snyder and Hughes. 

It is convenient to introduce a reference system, say 
K(xyz), associated with the symmetry properties of the crystal. 
For example, for NaNO, the z axis would be set along the 
@rystal triad“axis. It is’ also convenient to introduce a 
Second) reference system, Say Ke(x vy Zz ), withez, “coinciding 
with’ the magnetic field Hy direction. The polar angles of 
the z” axis with respect to the crystal reference frame K are 
denoted by @ and ¢. The quantities BY and Qh in equation 
(II 69) can be readily transformed from one reference frame 
to another with the use of a suitable rotation operator, the 
Euler angles for such a transformation being a = $¢, 8 = 6 and 
vy2=e6Q (Hdmonds, 2957) 7*°rollowing the. tsual procedure: (Van 


Kranendonk 1954, Pietila 1968) the quadrupole Hamiltonian 
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can be written as 


Hee ee Oy eps. (6, a Bo 
Ue ena 


ial 72) 


The term ols , an irreducible tensor, is in the K~ reference 
system representing the nuclear quadrupole term in a manner 
that is. independent of 8: and @. lLakewrse, the term BY is an 
irreducible tensor which is independent of the angular para- 
meters and is defined in the K reference system. The angular 
parameters are contained in the rotation operator which trans- 
forms BY into the K’ system, thereby putting the whole Hamil- 
tonian into that system. The matrix elements for the spin 


states are given by the expressions 


1 

ena <m+1|Q*+|m> = Hi /62 el ome.) Cl en een trem. 3) 
. 
2 aia <m#2|Q7*|m> = (76/2) k[ (Im) (I¢m-1) (I4m+1) (I+m+2) ]? 


(II 74) 


Considering a relaxation process due to quadrupole 
interactions (Pietila 1968, Van Kranendonk and Walker 1968, 
Snyder and Hughes 1970), the probability per unit time of a 
transition m, n-m+tu, n , where n and’n° are the initial and 
final lattice states whilst m and m+u are the initial and 


fined spin states, is given. to first order, by 
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W = 27 
m,n-m+u,n~ <a 


The 6 -function imposes energy conservation constraints 
on the transition. If the lattice is in thermal equalibriun, 
then the probability of spin transitions m>m+u can be re- 
written by replacing the phonon states by their thermal 
averages “This ;"theysprn “transition “probabrirty per Puntrt “time 


becomes in such a case 


Ee ils | Ora eae --D + DE-- Cu" uum). (It 176) 
where 
c(u",u>",u,m) = 27 5 <n?|BY )#[n>*<n-| (Be) nee 
R nn 
x 6(E. - + Pasian E.7 En’ (leo 73) 


Let us first suppose that the transition probabilities 
Wy and Wo depend on the initial spin state (Snyder and Hughes 
1970). Such a consideration was discussed earlier in connec- 
tion with the evaluation of the rate equations for the I = 3/2 


Spin system. We therefore rewrite Wi and W. as wr and We 


respectively. The two transition probabilities then become 
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Wesmay = Wy (2m+L)*(I4m+1) (I-m)/[2I(2I-1) 7] (II 78) 
and 
Wsmeg = Wo(I-m) (I-m-1) (I+m+1) (I+m+2)/[21(2I-1)7}. (II 79) 


If the transition probabilities do not depend on the 
m value, then the above equations reduce to the form of the 
expressions given by Yosida and Moriya (1956). Using equations 
(EL 43) POCLE 74) GandSCil 476), the spin transition probability 


can therefore be written in the form 
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By substituting for BY in (II 77) it can be shown (cf. 


Pietila 1968) that Yeh and Ww can be expressed as linear 


combinations of the quantities 
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where @,8,0 ,6 are the Cartesian, cordinates. x,y,z, which.for 
future simplicity we shall lJabel.as 1,2,;3. The quantity 
Mi g0787 is a component of the fourth rank=tensor representing 


the fluctuations in the electric field gradient. In the 
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general, this tensor will have 81 components. Using equation 
(II 80) and the symmetry properties of quadrupole interactions, 
it can be shown (Snyder and Hughes, 1970) that the M-tensor 


Le.reals Liwenen. rOldiows Lromy ii. 81). that 


and from Laplace's equation that 
Py Cee en = L- Mee = 0. lain o:3)) 


This reduces the number of independent tensor components 

from 81 to 15. The number of independent components may be 
further reduced depending on the particular crystal symmetry 
and this has been carried through for the 32 point groups by 
Snyder and Hughes (1970). Restrictions on the sign and mag- 
nitude of the M-components can also be obtained from equation 
(II 81) and it can be shown that these are (Snyder and Hughes, 


1970) 
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For the indirect harmonic Raman process and the 
anharmonic Raman process, the energy difference (Batu em? is 
small compared to the energy of the phonons and may be neg- 
Nected = inwequation (Li 961 ).).Lt “follows that: the M-components 
and hence the transition probabilities are independent of the 
parameters wiand m. In the case of the direct process the 
change in the nuclear spin energy is equal to the phonon 
energy and this approximation cannot be made. The M-components 
and ne are therefore dependent on u and m. 

Other relaxation mechanisms such as diffusion and 
torsional foscillatitonwof molecular.groups also fall into the 
former category where the M-components are independent of the 
parameters u and m, so long as the correlation time is suf- 
ficiently short. For those group motions where the correlation 
time is long, then the M-components will depend on uw and m 
(Snyder and Hughes, 1970). In practice, even where the M- 
components do depend on u and m, the m-dependence may be negli- 
gibly small_if the static quadrupole interaction is much less 
than the Zeeman interaction. 

Finally we-note, that: the orientation dependence: of Wy 
and W, is independent of the mechanism that produces the quad- 


rupole relaxation. For the case of a single crystal possess-— 


ing triad symmetry, the quadrupole relaxation transition prob- 
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abilities are given by Pietila, 1968, Snyder and Hughes, 1970) 
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SECTION III Apparatus and Experimentation 


a) Spectrometer 


The double resonance experiment that we wish to 
carry out requires power being applied simultaneously to 
two different frequency components of the nuclear magnetic 
resonance spectrum. This means therefore that the supply 
of power to the spin system and extraction of the infor- 
mation pertaining to the resonance must be achieved by a 
Single coil, leaving the second orthogonal coil to supply 
power to the second spectrum component. Crossed-coil 
spectrometers such as the Varian Associates models are 
therefore unsuitable for this type of experiment. Spec- 
trometers of the bridge network type are single coil units. 
However, they tend to be somewhat unstable and to respond 
to both the absorption and dispersive modes of the nuclear 
magnetic resonance signal, an undesirable complication. 

We initially constructed a Pound-Knight-Watkins single coil 
spectrometer (Watkins and Pound, 1951 and Knight, 1950). 
However, this unit could not be made to operate at the power 


levels that we required,with a sufficiently good signal-to- 
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noise ratio (Reed, 1967). ~The work reported in this thesis 
was carried out using a Robinson nuclear spectrometer, 
(Robinson, 1959) whose design was similar to that built by 
HOowLing (1966): 

In this type of spectrometer the nuclear spin 
system under investigation is placed in the coil of the tank 
circuit of an oscillator. The energy absorbed by the spin 
system during resonance originates in this unit so that 
Spin absorption will directly affect the performance of the 
oscillator. Therefore, by monitoring the oscillator level 
one can obtain information about the nuclear absorption of 
the spin species. 

The-tank carcuLt.of the oscillator. isha jparaliice! 
combination of inductive and capacitative elements, plus 
a resistive component that may be considered to be in series 
with the inductive part. The sample is placed in the in- 
ductive coil therefore exposing it to the circularly polar- 
ised components of the radio-frequency field resulting from 
the carrier developed by the oscillator. 

The absorption of energy by the spin system at 
resonance is equivalent to the introduction of an additional 


resastive. component, into the tank.circuit.. The resutting 
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Caption III a 


Figure III a shows the circuit diagram of the Robinson 


nuclear magnetic resonance spectrometer (Robinson, 
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echange*+in_ the*Q-factor of the 9tank-4ircuit-will thén*offer 
a means of monitoring the resonance. The Q-factor of such 


a parallel circuit is given by the equation 


(LILI ~-T) 


where L is the inductance, R is the resistance and Wo is the 
resonant pulsatance given by the condition w LC=1 where C 
is the capacitative element. This introduction of the addi- 
tional resistive term decreases the Q-factor,and the magni- 
tude’ of the. oscillation leéevel.falls. It is -the monitoring 
of this amplitude change that will allow the extraction of 
the information from the spin system. 

We shall now consider the particular design charac- 
teristics of the Robinson spectrometer whose circuit dia- 
gram.is shown in Figure III a . The spectrometer may be 
considered as consisting of two parts, a radio-frequency 
oscillator section and a demodulation-audioamplification 
section. 

Three 6688 tubes are employed in the oscillator 


section, which also performs the additional. role. of ampli- 


fying the radio-frequency oscillation by a factor of some 
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30.0 2.4 athe. fanal..stage of the. oscildator unit)is va: 6AK5 
mube eincorporated -in. ca limiter icincuit) . The feedback ais 
taken from the. output of the limiter and coupled onto the 
tank. Ci reuid: of the oscillator. The ALimiter outputs. a 
constant value that is independent of the input amplitude 
providing that this: input.is sufficiently large... Since.the 
preceding, section has a gain of 300,, this condition is,.easily 
satisfied. The feedback appears to the tank circuit as a 
constant current source with the degree of feedback being 
largely independent of the gain fluctuations in the early 
stages of the oscillator. While the characteristics of the 
limiter were adjusted to have the idealised square wave 
response, such a Situation was not achieved and a degree of 
ripple was always noticeable on the waveform. It is safe to 
say that the excellent performance of this spectrometer, 
particularly at the lower levels of power, is due to the in- 
clusion and successful operation of this limiter. 

Because of the signal-to-noise characteristics. of 


aS 


the Na NMR in NaNO one cannot observe the signal di- 


af 
rectly say by the use of an oscilloscope. Fast-passage 
techniques (Abragam and Proctor, 1958) are also unsuitable 


because of the frequency—dependence of the balance-mode 


system. Fast passage by variation of the magnetic field 
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is also ruled out since this would require the tracking of 
the saturation signal as the field is swept. We therefore 
employed phase-sensitive detection methods in conjunction 
with an audio-modulation of the steady magnetic field, the 
modulation amplitude being small compared to the linewidth. 
The result is an audio modulation of the radio-frequency 
carrier level of the order of 1 part in 10> Which is pro= 
portional to the first derivative of the absorption line. 

The audio-modulated radio-frequency carrier is 
taken from the last of the 6688 tubes of the oscillator sec- 
tion and is then fed into the demodulation-amplification 
unit. Here the signal is chopped, the radio-frequency compo- 
nents smoothed out by a series of low-pass filters, and 
the final audio signal is amplified by a factor of some 400. 
Thus the spectrometer output is an audio signal whose am- 
plitude contains the information about the nuclear spin 
absorption. 

Two other parts are. incorporated in the spectro- 
meter design although they are not essential to the obser- 
VatLon-ob the signal. One Ys a Circuit Ssultable for reed= 
ing the oscillator radio-frequency signal ‘into a frequency 
counter. The other is a calibrator unit (Watkins and 


Pound, 1951, Watkins, 1952). This unit G€mploys ad 6C4 "cube 
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Operated.in a nonlinear part.of its characteristic. The 
Principle here is) that, under such conditions, the plate 
resistance of the tube will be dependent on the instan- 
beous plate voltage... Therefore, if the, grid .veltage is 
varied at an audio frequency, we have introduced an audio= 
modulated resistance into the circuit. Furthermore, if 
the audio frequency is the same as that of the chosen 
audio field modulation frequency, then the detection and 
amplification units respond to it in the same way as they 
do to the results of the field modulation. The load of the 
tube is weakly coupled to the tank circuit so that we have 
introduced a pseudo resonance into the spectrometer. This 
circuit therefore offers a means of calibrating the sensi- 
tivity of the spectrometer's response to NMR and also a 
means of checking the linearity of the system as a whole. 
Care was taken in the design and construction of the 
spectrometer to reduce microphonics and pick-up, both of 
which.are of particular concern in. NMR spectrometers, par=- 
ticularly those operating at. low levels. The body of the 
spectrometer was made from 1/4-inch brass plate welded 
throughout. The body was divided into compartments each 
containing specific parts of the electronics of the spec- 


trometer thereby ensuring isolation of the various stages 
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Caption fire 


Figure III b shows typical resonances obtained during the 
main experiment. Those shown are for spectrometer observing 
powers corresponding to 20 millivolts peak-to-peak and 60 
millivolts peak-to-peak measured across the sample coil. For 


each level we show an enhanced and an unenhanced resonance. 
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from one another. The entire back plate of the spectro- 
meter could be removed to provide accessibility to all the 
Cancuitry. 

The. capacitance of the tank circuit, housed in a 
separate part of-the spectrometer away from the electronics, 
comprises a large variable capacitance of 50-600uuf and a 
smaller one of 5-25uyuf. The two are arranged in a parallel, 
the smaller one being driven by a synchronous motor via a 
set of reduction gears to provide the frequency sweep of the 
spectrometer. By using different gears, the sweep rate 
could be varied, a typical sweep rate corresponding 10kHz 
per: hour. Examples of typical resonances are shown in. Figure 


Tele wae 


b) The probe 


In order to. carry out. double-resonance. experiments 
it was necessary to construct a suitable crossed-coil probe. 
Such a probe must house the spectrometer and saturating 
coils as well as the field-modulation coils. The construc- 
tion of the probe should be such that microphonics are re- 
duced to a minimum and the components are well shielded. 


The probe was therefore milled from a solid block 
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of aluminium and all the coils were wound on teflon for- 

mers which fitted snugly into milled cavities of the probe 
body. The spectrometer coil was wound into grooves machined 
on a cylindrical former. The internal diameter of the cyl- 
inder and the length of the coil were chosen so as to provide 
a good filling factor for the particular size sample that we 
had avilable. The inductance of the coaxial cable that 
connects this coil to the capacitance part of the tank cir- 
cuit must be negligible compared to the inductance of the 
sample coil. There were a number of other considerations 
coneerning the sample coll? Firstly; ther Q-factor of* the 
tank circuit must be as large as possible in order to have 
maximum sensitivity in detection of the NMR signal. Secondly, 
the inductance had to be such that the resonance frequency 
was around 7.5 MHz, the 23Na resonance frequency in the mag- 
netic field that we had available. 

The spectrometer coil had an inductance of 2.1uH, a 
self-capacitance of 40uuf and a Q-factor of 135 at a frequency 
of 7.5 MHz when totally enclosed in the aluminium body of the 
probe. 

The single crystal of NaNO, was housed in a cavity 
in the end of a teflon rod which fitted snugly into the ver- 


tically-mounted spectrometer coil. The axis of the coil and 
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rod coincided, and»the»rod was free, to.rotate. about: the 
common axis. The top of the rod hada 360° angular scale 
with the associated vernier being mounted rigidly on the 

body of the probe. In this way, therefore, the orientation 
of the crystal with respect to the horizontal steady magnetic 
field, could.be)varied and, measured.to within Wiha: 

The saturating coil was wound in two slots milled 
into a teflon former. Between these two slots a hole was 
drilled to take the spectrometer coil. By inserting the 
Saturation coil into the probe in the horizontal plane, we 
could arrange for the steady field, the spectrometer radio- 
frequency field and the saturation field to be mutually 
orthogonal. The inductance of the saturating coil was 4.7uH. 
This allowed a sufficiently large saturation field to be 
produced at the correct frequency, with the coil being part 
of. the tuned. circuit of the saturation signal. amplifier unit. 

The modulation field was provided by a Helmholtz coil 
arrangement, the coils of which were mounted on the probe 
such that the field was parallel to the steady magnetic field 
direction. 

It is impossible to mount,the spectrometer and satu- 
ration coils:exactly orthogonal to each other to. avoid direct 


coupling between them. In order to overcome the coupling, it 
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was necessary to incorporate a balance-mode system the details 
©f_which.are discussed. in. Section III i) 6. It should be 
noted that the inclusion of this network will impose further 
constraints on the values of the inductances and Q-factors 
Orsthertwo coils: 

The probe was mounted on rails between the pole 
pieces of the magnet. This allowed the probe to be moved in 
the gap and scales attached to the arrangement allowed the 


fixing of the probes position in the field. 


c) Phase-sensitive detector 


The final output from the Robinson spectrometer con- 
sists of an audio-frequency signal whose periodicity is equal 
to that of the field modulation so long as the modulation 
amplitude is small compared to the linewidth. The amplitude 
of this audio Signal is then proportional to the first. deriv- 
ative ot the.nuclear absorption lineshape. If the spectro— 
meter is swept through the resonance, the spectrometer output 
gives. a representation of the first derivative of the line- 
shape (Andrew, 1955). The phase-sensitive detector is a 


device that produces a d.c. voltage proportional to Acose 
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where A is the amplitude of the Fourier component at the 
frequency of a predetermined reference signal and ¢€ is the 
phase difference between the component and the reference 
Signal. By feeding the reference signal and the modulation 
coils from the same-audio oscillator and applying the spec- 
Peonocer OULDUL tor the detector input, the dices output will 
be the. first derivative of the NMR signal. This output. can 
then be displayed on a strip=—chart recorder. 

Our phase-sensitive detector was built according to 
the circuit diagram given by Schuster (1951). Preceding the 
mixer stage was a twin-T narrow band amplifier tuned to the 
40HZ modulation frequency. This narrow band amplifier elim- 
inates unwanted harmonics of the modulation frequency at 80Hz, 
etc, which will otherwise also give a d.c. output. These 
other harmonics will be particularly present in the output 
of the spectrometer if the modulation amplitude is not small 
compared to the linewidth of the resonance. The narrow band 
amplifier also cuts down the bandwidth of the noise entering 
the detector. The bandwidth of the amplifier depends on the 
Q-factor of the twin-T network. The Q-factor was chosen to 
be approximately 10, since a much higher Q imposes too strong 
a. demand on the stability of the audio-oscillator frequency 


output. A low-pass filter with a time constant that could be 
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set at values from 1 to 80 seconds followed the Schuster 


mexer CLYCULt. 


dad) Saturating equipment 


The saturating Signal was obtained from a Marconi 
radio-frequency signal generator model TF2002. The frequency 
stability of this unit was a few parts in Mo per hour.) The 
frequency of the signal generator was monitored by a T.S.I. 
frequency counter, model’ number (“385R ;which had an accuracy of 
3 parts in chi A vernier frequency adjustment on the signal 
generator enabled the frequency to be set to within 10 Hz at 
Wey, MHZ 

The output. from the signal generator was. fed into a 
tuned radio-frequency amplifier developing a maximum voltage 
of 80 volts peak-to-peak. The output from this amplifier was 
fed into the saturating coil of the probe and a portion was 
also fed into the flux-balance system as discussed in Section 
ETT a). 

The Marconi generator had facilities for amplitude 
modulating the radio-frequency signal with known modulation 


depth. This facility was of great use in the initial inves 
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tigation.of the.functioning of the spectrometer. 


e) The Magnet 


The large magnetic. Field required to lift the nu- 
clear spin degeneracy was provided by a permanent magnet 
which had eight-inch pole pieces and a two-inch gap. It 
produced a field of 6700G and, as is usually the case with 
permanent magnets, its field was susceptible to thermal 
variations. This magnet had a temperature coefficient of 
approximately 1 gauss per °c. In order to overcome this 
problem, the magnet was completely covered with one-inch 
thick sheets of styrofoam which were glued to the yoke of 
the magnet. This arrangement was then covered with a wooden 
box which had a removable strip to allow access to the gap. 
This system therefore minimized the effect of variations in 
the temperature of the laboratory on the magnetic field 
strength. Furthermore, the laboratory itself was thermo- 
statically controlled and air-conditioned. 

The magnet was situated on a concrete slab laid in 
the floor of the laboratory. However, as we shall see in the 


discussion of the flux-balance system in Section III g), the 
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problem of vibrations originating outside the laboratory 
remained an inconvenience in the final experiment. 

An aluminium shelf was bolted across the magnet yoke 
near the magnet gap to support the flux-balance system. The 
Robinson spectrometer was placed on top of the magnet in order 
to minimize the length of the spectrometer-to-probe coaxial 
cable. However, care was taken to keep the spectrometer away 
from the fringe field of the magnet. 

The homogeneity of the magnetic field was investigated 
using an NMR sample comprising a saturated aqueous solution 
OteNaNO.. .One, tenth, of a.ce..of this: solution, was placed.in 


3 

geduaktZepiitad... Avlitite FeCl, was added. to: the NaNO, Sonu 
tion so that the Bei Lens would reduce the spin-lattice 
relaxation time Ty of the sample. This allows the use of 
stronger spectrometer observing power so as to improve the 
signal—ko—noise- ratio of the,.resonance. 

Detailed homogeneity plots were carried out in the 
gap of the magnet. Although the field was not symmetric in 

5 


the. gap,,aeregion with a, homegeneity of about 1 part. in 10 


over lcc was found. 
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f) Frequency Marking System 


The resonances were recorded by sweeping the spec- 
trometer frequency as opposed to sweeping the magnetic field 
which would require simultaneous frequency—sweeping of the 
Saturating Signal. Since we wanted to set the saturation 
frequency to within 30 Hz of the centre frequency of the 
Spectrum component being saturated, we needed to be able to 
measure the resonant frequency to this kind of accuracy. The 
spectrometer frequency sweep was not perfectly uniform. It 
was therefore not possible to use a frequency counter since 
this could not sample and display the varying frequency with 
enough accuracy without a substantial time lag. We now de- 
scribe in some detail the frequency=marking procedure adopted, 
since it is our impression that such a scheme is neither 
widely used nor easily understood. 

In essence, the system consisted of the following: 

A crystal=controlled RMS frequency standard, model number 
M76BR,delivering a signal at 100kHz which is accurate to 1 
Pati. rin toes This signal is then fed into a series of decade 
divider and multiplier. circuits and“amplifier units. producing 
Signats at iMuz, LOOKHZ, @0kKHz2, IkNz and L0UH2Z, DOEm in tne 


form of sinuSoidal fundamentals and pulses. The 10kHz pulse 
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was then combined with the 1MHz or 100kHz fundamental to pro- 
duce a series of Signals 10kHz apart throughout the radio- 
frequency band. This mixing and dividing was done using an 
Airmec frequency-standard unit. 

The spectrometer signal was picked up using a lead 
that was inserted into one of the oscillator tube cans on the 
spectrometers The spectrometer’ signal”-and the mixer output 
were then fed into an Eddystone communications receiver which 
had a band-pass of 5kHz. Tuning the receiver to the spectro- 
meter frequency meant that the audio output from the receiver 
consisted of the beat note between the spectrometer and the 
nearest 10kHz component marker. As the spectrometer sweeps 
through a range, the frequency of this beat note will change, 
decreasing as the spectrometer sweeps towards the pulse marker, 
and increasing as it sweeps away from the pulse marker. The 
beat note was fed into a loudspeaker and also to the X-plates 
of a double-beam oscilloscope. The 1kHz fundamental signal 
from the Airmec unit was fed to one set of the Y-plates and 
a *Vroon? Fundamental ’sronal was fed “to the-other -secr*or*y-= 
plates. ‘The resulting Lissajous figures indicate changes in 
the spectrometer frequency of 100Hz and 1kHz increments with 
respect to the 10kHz marker. After some practice, it was 


possible to identify these figures. Using an event marker 
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Prequency “markers could” be’ put. on’ the recorder chart. 

What remains now is to identify the marker pulses 
which will allow the exact determination of the spectrometer 
frequency. This was done by first feeding a 100kHz pulse-to 
the receiver. As»the receiver tuning control is varied, the 
Signal level meter (s-meter) on the receiver gives a large 
deflection each time the receiver is tuned to a 100kHz 
harmonic. Thus, having calibrated the receiver, it is pos- 
Sible to dial the receiver to within 10kHz of the desired 
frequency by simply reading the receiver scales since the band— 
pass of the instrument is only 5kHz. The tuning capacitances 
of the spectrometer are now adjusted to tune in the spectro- 
meter to the frequency set on the receiver. The 100kHz pulse 
fed to the receiver is now replaced by a 10kHz pulse. Thus, 
the spectrometer is now beating against a 10kHz pulse whose 
frequency is known,with the resultant Lissajous kHz pattern 
going 5,4,3,2,1 as the spectrometer sweeps Cowards ‘the marker 
and 1,2,3,4,5 as it passes the pulse and.sweeps away. The 
process would of course be reversed if the spectrometer was 
in the frequency decreasing mode. Thus, by observing the 
pattern, markers can be put onto the recorder trace every 
1000Hz and specifically identified. Observing the 100Hz 


Lissajous pattern will allow markers to be put on at 100 Hz 
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intervals. Finally, the 100Hz divisions can be subdivided 
by eye to give a final indentified frequency accurate to a 


few tens of Hz. 


g) The Flux-balance system 


In the performance of a double resonance experiment, 
meaningful results can only be obtained if the direct pick- 
up between the saturating and spectrometer signal is negli- 
gible. This is because the NMR sensitivity of the. spectro- 
meter depends. on the oscillation conditions of.the spectro- 
meter. and this can be changed by an external signal, particu- 
larly if the external signal frequency and the spectrometer 
frequency are close. 

In an idealized case of the saturating coil and the 
spectrometer coil being exactly orthogonal to one another, 
no pick-up between the coils would result. However, as was 
pointed out an the discussion of the probe, these coils cannot 
be mounted orthogonal to one another to the required accuracy. 
Furthermore, the coils are by necessity in a probe whose 
dimensions are restricted by the dimensions of the magnet gap. 
There will in general be a signal induced in the sample coil 


from. the saturating coil via. the body of the probe. This 
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Caption Tift c 


Figure III c is a schematic representation of the orthogonal 
and in-phase pick-up components considered in the evaluation 


of. the flux-balance system. 
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would in fact exist even if the two coils were exactly 
orthogonal, so long as the probe body is not symmetrically 
located with respect to the coils. 

It was pointed out in Section I that the saturating 
Signal had to be at least 40 volts peak-to-peak in order to 
ensure essentially complete saturation of a particular compo- 
nent of the spectrum. Also, the spectrometer level was in 
the range 10 to 100 millivolts peak-to-peak, where the en- 
hancement is a linear function of the observing power. It is 
obvious from the difference between these two magnitudes and 
the fact that one wishes to make measurements at crystal 
orientations were the two frequencies are as close as 40kHz, 
that a sensitive flux-balance system is required. 

Let us first. consider the form of the voltage pick-up 
in the spectrometer coil due to the signal in the saturating 
coil. The resultant pick-up may be considered as having an 
amplitude Corivas 1s allustrated ian Figure f2i ces “Ths 
pick-up can be resolved into two components one in-phase with 
the saturating voltage and a smaller component orthogonal to 
the saturating voltage. The in-phase component is due to 
direct coupling between the coils, and the orthogonal portion 
is due to induced pick-up via the. resistive body of the probe. 


The in-phase component is balanced out using the 


: cal ae 
er Siew efioo ows oth gs nove jax 008 , 


ne ylleaistomnys Zon at vbod side edt 6 enol Ce, 
v | | -2fbos ani od segues astm 
) ; PET aRe ont gent I ihaal at jy! besniog eaw a1 
7 ot tsbxo ft iseq-os -*tesg AS slov Oh tenet 46 ed oan 
i ,-oqmod ts hvoks %8g, 5 io noitsigtee ayeiqnoa vitshineeee & a0 
i at esw Lovet sodonorseege arid aka, MEASISGS silt 20 de 

a a 


Hens ois sxsw re ot-1BSq etlovilLim ial os OL 


3 FN 7 5 ‘ 
st i ‘ee ay P« 5 = 2 
b at tnom 


, 
| 


eb 3D). .zewog piiivreede eit’ to nokdona® : BPI 
_ bats | HR Nel arti ows 9aedt neawiad eonoxott ib: srt mor 

fsAyeyso ts adtasmesvesem sdem od eetelw ero seg) goats 
.Sht0h as se0l> 25 sis sotonsapers ows ods Si9W anotss: 


| | boxiwpen, ai meseye sons Ladhxiil® ovis iense a 

ae gu-wolg apetiov eft to mao? edd | 

Pima antec ors ne tenpte ang ey ‘oub fieo xedomossoege @ 

6 pi ived an poxsbienes od em: qu-xoig vasdiuaoe ant 

| ‘end, mo gs 0 e7pet ‘ei bdterdeulLi ait as : Lhe ebusts 
giiw seprg- ii end " ednsdogmos ows osme pevleass od. ano. worse 

‘i od. penoeers ts dnenogmos xolieme 8 Bis svediov | sa 


od sub at  silenigtinady seeda~at ont. ey pais 


15 


Caption: ELI id 


Figure III d is a schematic diagram of the flux-balance system. 
Particular attention is paid in the diagram to the grounding 


and shielding of the units. 
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mutual inductance coupling between a loop inserted in the 
ground return of each of the two coils as shown in Figure 
III d. These loops are single turns thereby keeping to a 
minimum the additional inductance that they will introduce 
into. their respective circuits. This-is of particular im- 
portance in the case of the spectrometer coil since it is 
imperative that most of the inductance in this line is lo- 
Gatedwat j thesecrystal . -Furthermoere; »theinelusion;jo£ the-Loop 
must. not markedly reduce the:Q-factor»of»the,tankscireuit. 
The two loops come within a few millimeters of each other and 
are supported along their whole length in the housing by 
styrofoam chips to minimize microphonics. The degree of 
coupling between the two loops is varied by means of a semi- 
circular vane that is rotated about an axis parallel to the 
common axis of the loops. This vane is made from one six- 
teenth brass plate and is one and a half inches long. By 
choosing the correct sense of the winding of these two loops, 
one will introduce a mutual inductance into the circuit that 
is opposite in sign to the mutual inductance between the spec- 
trometer and saturating coils. Furthermore, by altering the 
position of the vane between the two loops, the magnitude of 
the coupling can be varied until complete cancellation has 


been achieved. The position of the vane was varied by means 
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of a 2000:1 reduction gear, such a reduction being necessary 
to provide the required fine adjustment. The use of a mutual 
inductance coupling introduced into the ground—-returns of the 
two coils was first reported by Blume (1962). 

In his paper, Blume also suggested a means of re- 
moving the orthogonal component. Blume's suggestion is to 
feed the spectrometer coil with a signal of variable phase 
and amplitude which is coherent with the saturating signal. 
This involves applying a portion of the saturating signal 
through a phase—shift network to the spectrometer coil via a 
potential divider. In our particular experiment, we require 
that the tank circuit has a large Q-factor and that this will 
not be significantly reduced by the orthogonal balance-mode. 
Furthermore, aS was stated previously, the voltage applied to 
the saturating coil must be at least 40 volts peak-to-peak. 
These conditions could not be satisfied by the circuit sug- 
SestedcbytBtumen i Firstiyyointerdersthacsthe@-f£aceer is not 
significantly reduced, the orthogonal signal must be fed onto 
the spectrometer coil via a network whose effective resistance 
is much greater than the parallel resistance of the tank 
c@roeuwitewhi chil isvo€é- ches order of 10*ohms. Since the output 
impedance of the delay line is 450 ohms, a resistance of the 


order of 10> eohmsn iseéraéquiredeeithe result isithat theectmrent 
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reaching the spectrometer coil is insufficient to cancel out 
ene. orthogonal@picksiupl*” A°furthér: dbfFiculty arises at’ the 
input to the delay line. The source of the coherent signal 
fed into the orthogonal mode is the voltage amplifier that 
féeds (the saturating coil. «If connectéd directly to the delay- 
line, whose input impedance % 450 ohms, the saturating coil 
is overloaded by the orthogonal network with the result that 

a sufficiently large voltage cannot be developed in the sat- 
Unrating coil. -oThisidifficulty“could*only be avoided ‘by ‘in- 
cluding a resistance of the order of 10ohms in the delay line 
input circuit. This meant a further reduction in the magni- 
tude of the signal that finally reaches the spectrometer coil. 
Therefore, it was not possible to use this arrangement in 

our apparatus. 

The network we employed (Hughes and Reed, 1970) is 
schématically shown \in!Figure III-d-. The coherent ‘signal is 
fed from the amplifier into the delay line via a small 
Gapacitances*® The delay line is terminated by its character- 
istic impedance thereby providing a delay line output that is 
independent of the phase or delay line settings. The output 
is now coupled, via a 560 ohm resistance, into the tank circuit 
via an additional small loop introduced into the spectrometer 


coil ground-return lead as shown in Figure III d . The degree 


Me dup! e089 ot snsiociavend ak. Lio meveanransige! | 


ead 35 Roetrs yrivortE seit A iqbtoig Le 
i ‘dene te jnererdoo sit toss ‘eid: .gakl yelob oid 
| a  Geds roitiiqus aps low eds af siom Leropeds1o: 
jaye ioe eid’ oF Ylipextb betéeykoo 2D ites cnbsenusas ait 8 
a! ‘Lo paivatntce oft yamdo Ob8 0 sonskoumt suqetd «« 
| gene Sivesry ofa ddiw axnowten) hennaplatits eng yd Bb 

dee b4+ ab beqolevsh od tones: eens lov" eprei yiste 
Y i yYt hebiovs od vino olvtos ya Loo LTP: side 
cert velek octt mt endo Of Wa 1Shwo any to cbaddatees 5 | 
| barn sit ak nolitovbet asdza0 & FaRoR vain wtivont ay 
(hte xetomortosae ody eodgest yLieal® santa Lace be ong Bo 
: Five snonmaerAxts bked oan of sie tat ton: based 10% 


} re , * 
4 ~ Gi. ~ ae # 


ah 


ee (OTR! . beet bab asdguk) heyeLants uel faowion on - 
et Lanpie taowedoo sit | UID syupht ae ‘owed: y seme 
| flene « wiv sch yeleb oft cotati sous tens aca 
|) yMedtosxeds Bat vi Dedenlanes at omtl ysleb ea? «8 
|) Meg Barts tuqsu0 SAzL ysleb ® pabbiverg yours einciieiae 
7 | Aaggte ed? /epatssse daft ysis x6 semdq-eits 20 antiga 
Obiworto Anes edd oini ,sonsteiazox mdo- (082 5 iv belques mom tek 
Mi | gevenoxtosqe sft vin Beowhorsai qool LLsime Tsnoisibbe am Biv 
— .-& InT Saupls nt rt in NE. ie Coee 


79 


of coupling can be varied again by the use of a vane as in 
tChesecase of@theetin=pmaser couplingve Sdnceithesreturn! doophis 
only a single coil, the effect on the Q-factor and inductance 
of the spectrometer coil is minimal. We found that this 
arrangement allowed.us to cancel the orthogonal component 
without unduly loading the saturating coil. For the flux- 
balance system as a whole, we have independent control over 
the amplitude of the in-phase component, and independent 
control over both the phase and amplitude of the nominally 
orthogonal component. In principle, of course, one should be 
able to balance out the pick-up using only the orthogonal 
mode. However, the main component of the pick-up is the in- 
phase component caused by the mutual inductance coupling be- 
tween the spectrometer and saturating coils. A separate mode 
were this component is cancelled greatly reduces the demands 
on the phase-shift properties of the delay-line. It is 
therefore a great advantage to have these two independent 
balance modes. 

As is illustrated in Figure III d the orthogonal: and 
in-phase modes are housed in compartments separated by a 
perspex insulating sheet. Both compartments were rigidly 
constructed from one quarter inch brass plate to form a 


completely” closed unit.®%. This’ provided shielding’ and: reduced 
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the microphonic effects. The balance—mode system was located 
close to the spectrometer and probe to minimize the length 
of the connecting cables. 

The correct grounding and shielding of the apparatus 
zs or the greatest™-importance..and*this*is particularly true 
as regards the case of the probe and the balance modes. The 
Spectrometer and the saturating coils must have ground—return 
leads that are entirely separate so as to avoid interference 
between their respective signals. Furthermore, the housings 
should be grounded, yet in such a way that no ground—loop is 
created. The orthogonal-mode housing is grounded by the 
coaxial cable that brings in the signal from the saturation 
power amplifier. This casing also provides the ground for the 
delay line and in turn for the orthogonal coil that couples 
with the spectrometer coil. Thus the orthogonal—mode network 
is completely separate from the rest of: the circuitry, except 
for the inductive coupling, yet it is always enclosed ina 
ground casing. It-can be seen from Figure III-d) that two 
connections are necessary between the probe and the in-phase 
mode housing. To avoid the ground loop problem, and yet 
maintain complete shielding, the ground of one coaxial cable 
is not connected to the housing even though it completes the 


return path of the coils. This was achieved by mounting the 
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Caption III e 


Figure III e is a block diagram showing the interconnection 
of the units that-make up our apparatus. Each unit is 

labelled to signify whether or not it is grounded. Also the 
distribution rack feeding the particular piece of equipment 


is ishown. 
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coaxial connectors on insulating boards. This complicated 
ertuacion becomes clear by. consulting, Figure 1lilsdi.) The 
two spectrometer—coil leads from the probe are especially 
sensitive to stray radiation and these were doubly-shielded 
using braid slipped over the cables and grounded only at the 
probe end of the cables. Thus, we have two completely in- 
dependent ground—return paths for the two modes of the flux- 
balance system, and both are always contained in a shielded 
environment. 

The performance of the balance mode system is discussed 


Mieaetaid. on rSectionviliIl 1) 6. 


h) Interconnection of units 


We have already discussed in detail the most important 
parts of the equipment. During this discussion it was pointed 
out, particularly in connection with the balance-modes and 
the spectrometer, that grounding of the apparatus is of the 
utmost importance.. As can be seen from Figure III e¢ there 
are several. units, woich if frounded, would create: qrovnd— 
loops. 

As is shown in the figure, the main power lines come 


to the apparatus by way of- three distribution racks each with 
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Six outlets. The line voltage feeding particularly sensitive 
preces Of equipment is stabilized by a Tine régulator:; “We 
HRave-inctuded. in the figure an indication as to the distri- 
bution rack that feeds a particular unit, and whether or not 
the piece of apparatus is made floating with respect to the 


electrical ground, 


i) Preliminary experimental procedures 


This section is concerned with various measurements 
and calibrations that were necessary before the main experi- 


ment could. be carried out. 
1) The spectrometer 


The performance of the spectrometer as a radio-fre- 
quency amplifier and detector was first checked at 7.5MHz 
using an audiomodulated radiofrequency carrier from the 
Marconi signal generator. This signal was traced through the 
various sections of the spectrometer, and components were 
altered until the performance of the system proved satis- 
factory. ._ Further adjustments of the. spectrometer. components 


; : Za 
were carried out using the calibrator unit and the Na NMR 
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Signals from an aqueous solution of NaNO, and a single 


3 
Cpvstalsor NaNO.. Suffice to say that this procedure took a 


good deal of time before satisfactory resonances were obtained. 
2) Linearity and long-term stability of the spectrometer 


The linearity and long-term stability of the spectro- 
meter was investigated using the calibrator unit. A known 
calibrator voltage, sufficient to give a deflection on the 


chart comparable with that of the <3 ga Signal in NaNO was 


37 
fed into the spectrometer. The output from the phase-sensi- 
tive detector.was then recorded for several hours. The long- 
term stability of the deflection caused by the calibrator 
Signal was found to be satisfactory. The phase of the cali- 
brator Signal was then changed by 130s. resulting in a de- 
flection in the opposite sense on the chart recorder. The 
total deflection corresponding to the magnitude of the cali- 
brator signal was then measured. As pointed out in Section 

TII a), the calibrator unit: provides a pseudo resonance. The 
linearity of the spectrometer to NMR signals of different 
strength may therefore be checked by measuring the spectrometer 
sensitivity. to different.values of.the calibrator voltage. 


The procedure described above was therefore repeated 


for different calibrator voltages. The sensitivity in milli- 
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meters .of recorder deflectton*per*mPillivolt.of calibrator 
fignel 18 shown in) Table Illia, Lters apparent *that the 
system 1s-linear in:its sensitivity to-within the limits 


imposed by the signal-to-noise ratio. 
3) Phase-shift network 


The phase-shift network provides the coherent 40Hz 
Signals for the modulation coils, calibrator unit and the 
reference input to the phase-sensitive detector. Each output 
has an independent phase and amplitude control. Using an 
NMR sample containing an aqueous solution of NaNO., the phase 
of the reference signal was adjusted until a maximum resonance 
Signal was recorded. Then, using a calibrator signal, the 
phase of the calibrator was adjusted for maximum responce on 
the chart. <In«this: way, the¢correct iphase relationships be= 


tween the various signals was obtained. 
4) Field modulation 


The 40Hz modulation of the steady magnetic field is 
provided by means of the Helmholtz coil arrangement that is 
built dmto the. probe. A small. search-cotl was: used to4in= 
vestigate the magnitude and homogeneity of the modulation 


field produced in the region of the sample. The modulation 
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Table’ JiI.ea 


Calibrator signal Delfection Deflection/Calibrator 
Signal 
(millivolts) (centimeters) (centimeters/millivolts) 

~2 
96.0 4.64 483% 10 

134.3 per tak 2 77x 10 

-2 
174.0 8.43 4.84 X 10 
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field was calibrated by measuring the voltage induced in the 
search coil by a known modulation signal. The amplitude Me 
in volts of the induced voltage is related to the modulation 


amplitude in gauss by the equation 
Bea ay 
O ee 


NA W 
m 


Cia) 


where Wa is the modulation angular frequency, and NA is the 

area-turns of the search coil. Using an accurate flux-meter 
and a known magnetic field, the area-turns of our search—coil 
was found to be 103.8 cms*. The modulation field was found 

to be 1.6+ 0.1 gauss peak-to-peak for a 1.0 volt peak-to-peak 
output from the power amplifier that fed the Helmholtz coils. 
This measurement was carried out for various positions in the 
probe sample region but no detectable variation in the field 


could be found over a region of several ccs. 
5) Performance of the flux-balance system 


We have discussed in Section III g) the problems in- 
volved with the crossed-coil probe and we outlined the ar- 
rangement that was developed to counteract the problems in 
our case. In this section we shall be concerned with the 


performance of the balance modes in a double resonance exper- 
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iment. 

A gross indication of the saturating signal directly 
affecting the spectrometer can be seen by observing the spec- 
trometer radio-frequency level meter as the saturating signal 
is increased. The pick-up by the spectrometer coil results 
in the spectrometer being forced to oscillate at the satur- 
ating signal frequency rather than at its own free-running 
vawwe, ihe derector level indicates this by an initial 
Wectease due CO Suippress.on Of the free-runhang’ osc11 lation. 
followed by a rapid increase due to the complete breakdown of 
the spectrometer's operation. 

A more sensitive indication of the effect of the 
saturating signal on the spectrometer can be obtained by ob- 
serving the audio output from the spectrometer. This output 
contains a sinusodal component whose frequency is the differ- 
ence between the frequencies of the saturating and spectro- 
meter signals. As the influence of the saturating signal on 
the spectrometer increases, so does the amplitude of this 
audio output. 

The most sensitive criterion is the frequency of the 
spectrometer, particularly as illustrated by the Lissajous 
pattern on an oscilloscope. An extremely small unbalanced 


saturating signal. will cause a change in the Lissajous display, 
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without it being perceptible to the unmusical ear. 

The Mmost,important. criterion: is. the’ effect that the 
Saturating signal has on the spectrometer's sensitivity to 
an NMR signal. When the spectrometer is being to some degree 
driven, it-is found that the. NMR signal intensity is artifi- 
cially increased. 

This effect was investigated by measuring the inten- 
Sity of the enhanced signal as a function of the deliberately— 
intorduced coupling between the saturating and spectrometer 
coils. The crystal was orientated with respect to the field 
so that 9 = 45°. This corresponded to a frequency separation 
between the saturating and spectrometer signals of approxi- 
mately 42 kHz. A saturating voltage of 40 volts peak-to-peak 
was applied to the centre line component and observations of 
the high-frequency satellite intensity were made at an observ- 
ing power of 25 millivolts and.then repeated, for an observing 
power of 60 millivolts. The misbalance was measured in terms 
of the amplitude of the audio signal from the spectrometer 
output. This amplitude was measured in centimeters on an 
oscilloscope with a 25 volts/cm sensitivity, and.is to this 
extent an arbitrary measure. The date is given in Table 
Tti- be eand ashown, in Figure. li th.) For the casevorva, 2Zommi1 1 — 
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Table III b 


1) 25 millivolt level 


Misbalance 
(arbitrary units) (centimeters) 
a.O5 0.80 
J «25 E00 
7.07 ie. 510 
Lo LF 2.00 
ERPS) 30:0 
7.20 4.00 
7.40 6.00 
(Eee 7200 
Ta50 10.00 
Sie oD £3.00 
9.40 ESO 


9. 80 20.00 
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(erasamitias) 
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60 maLidavodte level 


(arbitrary units) 


Misbalance 


(centimeters 
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Caption?’ ILI *t 


Figure III £ shows the intensity of a satellite, 
recorded whilst the centre line component is severly 
Saturated, as a function of the deliberate misbalance flux- 
balance system. The intensity of the satellite is given in 
arbitrary units and the deliberate misbalance is expressed 
as centimeters of audio signal from the spectrometer audio 
output as discussed in the text. The variation of intensity 
is given for two values of the spectrometer observing power 
corresponding to 25 millivolts and 60 millivolts peak-to-peak 
acress the sample coil. For beth sets .of\data, the saturating 


and spectrometer frequencies differed by 40 kHz. 
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signal does not begin to increase until the pick-up ampli- 
tude exceeds 2cm on the oscilloscope. A change in the note 
does not become audible until the pick-up exceeds about 3cm. 
After 4cm, the intensity rapidly increases, clearly illus- 
trating that data taken above such a misbalance situation 
would lead to very erroneously values of the enhancement. 

It is noticed that the indication of the misbalance concerned 
with the suppression of the detector level and its subsequent 
increase is not observed until the degree of misbalance ex- 
ceeds some 17cm” Clearly at hiSucri terion TOotebalLance iis 
completely wrong and explains the large enhancement values re- 
ported in the author's M.Sc. thesis. For an observing power 
of 60 millivolts, the increase in the enhanced intensity 
occurs at a higher misbalance, as one would expect. 

For all the enhancement data reported in this thesis, 
the pick-up was maintained below 2cm on the oscilloscope. 

If, during, any run, the, pick-up exceeded this, value, the 
resonance was discarded and the run was repeated. 

Let us now consider in absolute terms the ability of 
the balance modes to cancel out pick-up between.the two coils. 
In order to measure the signal induced. at the saturating 
frequency in the spectrometer coil, it was necessary to use a 


communications receiver in order to reject 2nd and higher 


* We are grateful to M. Smith for pointing out that the 


2 
intensity « (pickup) 


ge ee Lay a 
~ i Leis, queisiy soft Pesos: sesszont os niped jon 
Sion edt at ‘spnisda’ A: _sqonset bigieo act bs fe) OR. | 
tot tuods ABSSOx5 gage Akg coy ba Liste sidibus. puaros: 
-agfit yltasio .asess omnd yibhast whisastrt, sd oh 
hottangie sonsisda ia 5 ove eveds: asts: Pr saad § 
.taomeonssine ait to sonisy vievontionss agtev od a6: 1 , 


wan 


benteonod So: ta lecdain edd to noLsegiead itt ad beoks 
Senge heise etx Bas Leva Let aIeh efts x noiessagqua ws 
re sansledeim to ss1pep ait Lisgan bevisedo Jon eh: t 

| af sonslgd) 102 nokveditto: e#as hase. ‘mort ¢ cook a 

“9x aeulav jnemeonsins spisl att. anislaxe: ‘bas ‘prow so 
tswog paivisedo iis tof .eiesds 98 we é* xottous ont ek | 
yo henodnt psonsdne ons at sanoioat: EF ha sebovtl 
stoscxs blvow sho S& sts. Ledi dn, sorigi = ga 

weherss eids nt hetzoqgss sjeb. Sidmeomade ont bia Brak a i 
ieqooeorLt foeo0.45da. ao mod wold oped idsindsit esw qu | 

capi sn tisy ais bsbesoxe qu-sotg ond fire. yas de 
sDetssqet asw gut eis bis bobuspeib eew 29 ; nt 

Bo yilids of? emiot ojtfords ak vebieso> won ew don 
,akion wd e8d asewiod quad tuo Lsotgs, oF eebom sousisd. 
(Pabdswwiss Sit 35 bestibet tsetse ort eusesom of 

| i sen 03) Yaseee0sa 2ew df fico 19345moxt0eqe ody mk 1 
eee ne bak joetex od tobvo at sevieosr amoitse 


alte 45: do pagtnteg 10% aoe -M od being sae all . 


of "| Sabaiaigy > > ppbaneeict 
=e AZ tt) a ney. 


94 


harmonics which otherwise would overwhelm the required signal. 
Because of the low input impedance of the receiver, a pre- 
amplifier unit was inserted between the coil and the receiver. 
This preamplifier had a low noise-figure, a gain of approxi- 
mately 200 at 7.5 MHz, and an input impedance of eo, ohms 
shunted by a 2uyuf. The preamplifier therefore did not 
Significantly load the spectrometer coil. 

A signal of 50 volts peak-to-peak was fed to the 
saturating coil with the flux-balance system disconnected, 
and the pick-up at the spectrometer coil was found to be 0.35 
volts peak-to-peak. (This signal was so large that it could 
be measured directly on a calibrated oscilloscope. This 
voltage is the sum of the direct pick-up due to the non- 
orthogonality. of the coils; sand the orthogonal component. 
caused by currents induced in the resistive body of the probe. 
The coupling between the saturating and spectrometer coils 
was then minimized using the in-phase flux-balance mode alone, 
and the pick-up was found to be 3.3 millivolts. Using both 
modes the minimum was 5 microvolts. The technique used to 
minimize the pick-up was to first use the in-phase mode and 
then make the final adjustment with the orthogonal mode. The 
in-phase mode was by far the most effective as one would ex- 


pect since most of the coupling is due to the direct coupling. 
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The frequency of the saturating signal was now 
changed over a range of approximately 5kHz and the resultant 
pick-up was measured by means of the preamplifier-communica- 
tions. receiver. combination.. By readjusting to a position of 
minimum pick-up it was found that the in-phase mode was 
frequency—dependent while the orthogonal was apparently fre- 
quency=independent. We attribute this to the effect of stray 
capacitances in the probe and the in-phase mode housing. The 
results of the frequency dependence studies are given in 
Hable ©li-~c and illustrated in Figure Lil-g- 

Let us now consider the temperature dependence of the 
flux-balance system. This was investigated by minimizing the 
pick-up and then changing the temperature of the laboratory 
by 5°C. Then the temperature of the laboratory was allowed to 
return to its original value, a process that took a few hours. 
Every ten to fifteen minutes, the degree of misbalance was 
measured using the centimeter criterion mentioned in Section 
III i) 5. At:-each temperature setting, the change in the 
saturation frequency needed to minimize the pick-up was found. 
Using the curve given in Figure III g, the frequency shifts 
were related to microvolt, levels and a temperature coeffi- 
cent of 20uvV/°Cc. was calculated. The data is given in Table III d. 


Let us now propose a model to analyse the frequency 
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Caption itd 


Figure III g shows the square of the pick-up voltage 
as a function of the saturating frequency as the saturating 
frequency is mis-set. The setting of the flux-balance system 


was unaltered during these runs. 
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Table III d 
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dependence of the pick-up of the system. Consider that there 
is a residual pick-up of. magnitude Vi that has a phase angle 
@ with respect to the in-phase component, and that this is 
frequency-independent. We also have an in-phase component of 
amplitude V>5 due to the direct coupling which the previous 
data shows is frequency dependent term. The resultant pick- 


up V is therefore related to these two terms by the express-— 


ion 
a oe 2 2 
[Vie] = (V+V,cos®) + (V)sin a) s GieLT 4.) 
Now we can write V5 in the approximate form 
V5 =K (v- we, (ClTee5)) 


where K is some constant, V is the saturating frequency and 
ae is the frequency at which the in-phase component V, 


vanishes, Thus we can write (III 4) as 


2 
2) = K2V7+ (2KV, cosd -2KV.) Vv + (K°VO-2K 
O O 


2 
‘ 5 al ab 
cosoV, Vy + Vi) ( 6) 


Such: a model would predict that the frequency dependence. of 
the pick-up would be in the form of a parabola, and this is 


indeed the case for the data shown in Figure III g). 
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Using a least squares analysis of the data, we find 


ne? eas (ara2 nh3s eyes! e(=1 sedoumeign O97 z6nesen(rTt 7) 


giving K the approximate value 35uV per kHz. The theoret- 
ical fit to the experimental data is shown in Figure III g. 
The question remains as to the causes of the pick-up 
Vi° By studying the operation of the spectrometer when 
recording enhanced signals it becomes apparent that the limi- 
tation of the flux-balance. performance is due to micro- 


phonics. We therefore atrribute Vi to microphonic effects. 
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6) Effect of mis-setting of saturating frequency 


To’ ensure the complete saturation of a component of 
the spectrum, the applied power must have sufficient magnitude 
and be applied at the correct frequency, corresponding to the 
maximum of the absorption curve. In the case of the recorded 
first derivitive of the resonance this frequency will be the 
crossover value. An insufficient magnitude or incorrect 
frequency results in the recorded enhancement being less than 
the value for a completely saturated component. The effect of 
insufficient power will be discussed later. In this section 
we shall discuss the effect of deliberately mis-setting the 
frequency of the saturating signal. 

According co. the BPP theory, the "saturation behaviour 
of a homogeneously—broadened resonance line is given by 
equation (II 65). However,it was pointed out, by Redfield 
(1955) that the BPP theory was invalid for S>>l, where 


2 


1 2 
07 a ieee, 2 


except for the particular case where the saturating power is 


and T. is the so-called spin-spin relaxation time, 
applied to the centre of the resonance. This was confirmed 
experimentally by Goldburg (1961) by measuring Mo the 

2 : 
z-component of the nuclear magnetization, of the 3Na spin 


system in NaCl as a function.of the mis-setting of. the satu- 
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rating frequency from the centre value. For S.= 1, the 
measured M, values were in agreement with the BPP theory. 
However, for S = 10, the BPP theory was obviously in: error 
except when the saturation frequency corresponded to the 
centre of the resonance, but the measured values were in a- 
G@reement with the Redfield theory. , According to this theory 
it is not possible to obtain complete saturation of a resonance 
(M=0) using a steady-state technique if the saturating fre- 
quency is mis-set. It was therefore vitally important to 
investigate the effect, on the satellite enhancement, of 
deliberately mis-setting the frequency of the saturating sig- 
Ralmappiied to thevcentre’ line. Our dateagis shown 1n- Figure 
Pie tOr a. field modulation amplitude of0..65 GCG and 1:30 (G. 
It is clear that mis-setting the saturating frequency seems to 
have a relatively small effect. Furthermore, the degree of 
mis-set necessary to produce a significant effect on the sat- 
ellite enhancement is dependent on the modulation amplitude. 
This effect can be explained by realizing that the field 
modulation has the effect of modulating the saturating frequency 
across the centre line. It thus appears that so long as the 
saturating frequency is swept by the modulation across the 
centre of the resonance then essentially complete saturation 


can be produced even though the saturating signal is mis-set. 
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Table III e 


i) Modulation amplitude 0.65G (0.73kHz) 


év Unsaturated intensity Saturated intensity Enhancement 
(Hz) (centimeters) (centimeters) 
0 6.600 Get 75 1.482 
= FOO 6.575 94525 1.450 
=250 6.875 10.200 1.483 
=5.00 6.575 95625 1.465 
= 150 6.775 9.600 1.418 
=1000 6.700 J+eL50 1.365 
=2250 6.700 8.750 e305 
— 2500 6.600 7«850 Lak90 
=h750 62825 Ws525 Le hOD 
-2000 6.675 7.000 1.048 
100 6.900 10.295 be AlS 
250 6.675 10.050 he DOZ 
500 6.670 9.900 1.478 
750 6.878 9.875 1.435 
1000 6.800 See) i e3Ge 
1250 7.000 921325 ARE IER A 
1500 6.725 8.275 L225 
1750 6.200) 1.100 Pelee 
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ii) Modulation amplitude 1.3G (1.46kHz) 


év Unsaturated intensity Saturated intensity Enhancement 
(Hz) (centimeters) (centimeters) 

0 6:10:25 95075 UA S10 
e530 5;800 Se175 i, 510 
7070 Spee: 8.875 1.498 
-1000 5.390 S26 1.493 
E1330 5. 90 8.275 Deeb 
=1670 Bip pe) Ae) WIENS) (raise 5 Yls 
=2000 5. 890 Te LO eee 

33.0 be 50 82 925 1499 

670 3 00 Bi750 1.483 

1000 62,050 9.000 1.487 

1330 Saree 8.300 Lai) 

1670 550 18 O12 1.324 


2000 64.075 Leone 6238 
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Caption PLl 


Figure III h shows the enhancement of the high fre- 
quency satellite due to severe saturation of the centre line 
as a function of the deliberate mis-set of the frequency 
of the saturating signal. The two curves correspond to 
field modulation amplitudes of 0.65 G and 1.30 G. The 
solid lines do not refer to any particular theory but are 


simply drawn through the points. 
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It should be possible to treat this effect by analyzing the 
frequency—modulated saturating signal in terms of its compo- 
nent sidebands. However, the important feature of the data 

is that it is not essential to maintain the saturating frequency 
at the centre frequency. In our determination ene en- 
hancement,the saturating frequency was always maintained with- 
2 30H, of the. centre of, the. resonance... Jt might, appear from 
Figure III h that the satellite enhancement is slightly de- 
pendent on the modulation amplitude at zero mis-set. However, 
later more careful measurements showed that this was not the 
case. 

The experimental procedure for determining the satu- 
rating frequency was to measure the cross-over frequency of 
the resonance to be saturated,in the frequency-—increasing and 
then frequency—decreasing sweep modes. These values were not 
exactly coincident because of the time constants of the appa- 
ratus. ft 1s Obvious that the correct cross-over frequency 
at which the saturating power should be applied is the mean 
of these two results. 

In order to check that no large errors were made in 
the determining of the saturation frequency, the enhancement 
of a component was measured for the case of the saturating 


power being deliberately mis-set by +lkHz. An incorrect 


ROL me 


, ext enisylsas. ed! soatte abide: saond oa stdieeog ed & 


-ogmed att to) emo? read Langie praventitee, Bepats 
sisbh ets Io oxvtbSst dnstrodHT oda a 


Yoroupet pnitsivise ods rube token, ae Ieitasess gon ai ‘em 
~f5 edt to robsamimredsb, ua fiz . yorsypedt nal 
 echtiw bentsiniem eyswis eb voneuipaad: gatietbice orld 3 
\Mexr2 ts0eqqs unpim 31 .soganoder ed Be patie aa to, 

| Lin old ribé Le at Jnemeonsdis etiliestsae oct ner (aa s 

| (sevewol ~tee-eim ones 36 ‘shagilome nbidesssbom, 8 Ht 


arty ton asw cind tsds bowode cane Sse igherso, 


rap 
v 


auden ois piinkmretsh xot aEnop era Lastniaeane | 
Bo yansupest? wSsvo-aa0xs: shit oundsem or sett 
Pn, bes eae yonsupet? sid at ebodssse ed’ een 
ton Stew esulsy eee |. s9bon ‘qsowe ontesexc 


ey re 


| -sqqs ect 2o adtnsdanos sms ont to seusobd 4 Biorntoc yitoss 
bie NVHeupos2 TSVO-B20T9 | Soanon edt tect stoi ai | salt ay a i 
_-. teem of? ef betiqgs sd, bivode aewog pobsnudne ait | 
ey astm owe | 
ai on ebem Siew axcxis eyzsl on 4ens toetg o xehs0 ae 
7 Saemeonadris ens  Yousupext noitsrutsa edd 20s | 
| paktexc2se ett to sesn odd x0} boaasem eew 


sGexr09e1 eh shale yd toe-eim erences ms es 


107 


saturation frequency would show up by the enhancements being 


different in the two cases. 
7) Calibration of the sample holder angular scale 


It was necessary to calibrate the angular reading 6° 
on the rotation scale fixed to the sample holder in terms 
of the angle 6 between the triad axis and the steady magnetic 
field. If the symmetry axis of the crystal is perpendicular 
to the. axis. of rotation, then 6 = 6° + y where 7 is a constant. 
The separation of the satellites is related to the angle 6 


by the. equation 


2 


ued, 3cos* eat) (III 8) 


2h 
where eQq/h is the quadrupole coupling constant. By meas- 
uring the satellite splitting as a function of 9°, we can fit 
the data to equation (III 8) to obtain the relation between 
6 and 6°. Furthermore, by comparing the measured quadrupole 
coupling constant with previous measured values (Pound, 1950, 
Andrew et al., 1962) we are able to check that the triad 
symmetry axis is essentially perpendicular to the axis of 
rotation ofthe. crystal. 


The satellite separation was measured at 30° intervals 
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Ghevatthrough asrangetof 180°. Such a procedure was carried 
out twice, once each for the two values of the azimuthal angle 
>. The point of repeating the measurements for the second 
setting was to check that the symmetry axis was still perpen- 
dicular to the rotation axis after the re-orientation of the 
sample to change the ¢ value. The data for the two cases is 
given in Table III f. Using a least squares analysis tech- 


nique (Whittaker and Robinson, 1925) a linearized equation 


giving the satellite separation £(0°) in kHzas. 


£(05) = 84.87+(-212..43)1cos26~+ (132.92) sin20~ (ITHT9) 


O 


can be obtained for the setting > = 2.81+1.5 For the case 


fe) : : ; 
the linearized form is 


Ge 3:26:02 150 
£0 ) ="84.07 (248.58) cos26°4(47.25) sin20, (ETE a0) 
The general form of this equation can be written as 


£(6°) = A + Bcos20~ + Csin26~ (Tr tale t) 


and the relationship between the angles 8 and 6° can be 


written as 


20° + 2y = 20 or Os 0° + 7. (TTL 2) 
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Therefore, we may write 


Rcos (20°+2y) = Rcos26°cos2y - Rsin2@“sin2y = Rcos26 Chir a 3) 


2 


where R = /(B? eeu es) 


If we write 
£(8°) = A + Rcoos(26~ + 2y) (Jit 14) 


then by comparing equations (III 11) and (III 13) we obtain 


the relationships 
Rcos2y = B, Rsin2y = -C and therefore Tan2y = -C/B (PETALS) 


By re-arranging equation (III 13), we can write 


£(9°) = A+Rcos(20°+2y) = A+Rcos20 = (A-R)+2R/3(3cos“6). (III-16) 
Defining X = A-R and Y = 2R/3 we have 
£(8) = Y (3cos*6 - X), emia aly 
. yA 


Knowing the values of A, B and C from the linearized equations 
we can calculate the angle 2y and the coefficients X and Y. 
O Om 
For the case ¢ = 2.8 + 1.5 we have y = 6 1° there- 


fore giving the relationship 


Bere roe Eee 1 (III 18) 
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Thus the satellite separation is given by 
Av = 167.1 (3cos76 — iM spsie wes sii cwas aby) Clepiar 19) 


The experimental error in the quantitles 165.6 and 167.1 is 
estimated at 1kHz. Equation (III-9) is thus of.the form 

Av a (3cos 0 -1). The quadrupole coupling constant is there- 
fore 334.0 + 2kHz. and is in good agreement with the values 
given by Pound (1950) and Andrew et al. (1962). 


O 


Por the Case’ o =)32.8 2.1.5 , by a similar technique 


we get 

@= 0° + 0° 37 (III 20) 
and a satellite separation given by 

Pee 6en0) (2cqs-6) — 169.0468. 0) (tao) 


In this case the quadrupole coupling constant is 336.0 + 2kHz, 
again in agreement with the accepted. value. We therefore 
conclude that the symmetry axis of the crystal was in both 
cases essentially perpendicular to the axis of rotation of 


the sample. 
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Captione Lil <i 


Figure III i shows the rhombohedral unit cell of 
sodium nitrate. For clarity only two groups of oxygen 


atoms are shown. 
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8) Crystallographic information on NaNO, 


Sodium nitrate crystallizes in the rhombohedral 
system (space group R3c, DS 4) with two molecules in the unit 
cell. As shown in Figure; lil,i;.there,are nitrogen atoms 
at each corner of the unit cell and also at the centre; 
the sodium atoms lie half-way between the nitrogen atoms on 


Ehnewtriad.axis.of.the,cell:at GNC Sao S LHe PDI Leaite 


444 
groups are planar, each nitrogen being situated at the 


dil 
44474 


centroid of an equilateral triangle of oxygen atoms whose 
place is normal to the triad axis. The sides of the oxygen 
triangles lie parallel to the projection of the cell edges 

on a plane normal to the triad axis, the orientation of the 
triangles alternating by 180° as one proceeds along the triad 


axis. The oxygen positions are given by the parameters. 


cad 
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The sodium atoms therefore have an environment of trigonal 
symmetry, and the axial field gradient tensor is axially 
symmetric. 

From close examination of the unit cell, it is clear 
that the environment of each sodium nucleus is similar. How- 


ever, the environment of adjacent sodium sites situated on 
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the. triad axis (as schematically shown in. Figure III 4) differ 
in one important respect. Since. the orientation of. the nitrate 
groups alternate by 180° as one moves along the triad axis 

the environments of the sodium sites are not identical. In 
fact each of the two sodium sites in the unit cell are in- 
equivalent, in that they differ by a rotation as well as a 
translation. In other words, the principal axis defining the 
symmetry of the sodium environments have common z axis along 
the triad axis but the x and y axis are differently orientated. 
This difference stems from the orientation of the oxygens 
associated with the nitrate groups. If the nitrate groups 
possessed azimuthal symmetry, as would be the case if the 
groups rotated rapidly about their own triad axes, the sodium 
Sites are no longer different. We note however, that since 


2 2 


a*vV/az- is the same at each sodium site and 3° V/0x™ = say ay 


ieev ace from Laplace's equation, the static quadrupole 
Brg ie is identical. for each 23a Site. This is why 3 
resonances and not 5 (2 pairs of satellites and a centre line) 
are observed. 

The question now arises as to whether there are two 
non-equivalent sodium sites in so far as quadrupole relaxation 


is concerned. The factor determining the quadrupole relaxation 


is the orientation of the externally applied magnetic field 
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Figure III j is a schematic representation of the 
two Na Sites believed to exist in sodium nitrate. The 
figure shows a side view illustrating the planes of nearest 
neighbour oxygen atoms that are associated with the two 23Na 
nuclear sites. Also plan views are given for the two sites 
showing the orientation of the nearest neighbour oxygen 
atoms with respect to the projections of the rhombohedral 


unit cell edges on the planes perpendicular to the triad axis 


Or che unit Cell. 
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relative to the environment of the relaxing nucleus. To 
consider this point we refer to the conventional xyz co- 
ordinate system for a rhombohedral system where the z axis 
coincides with the triad axis. and the x axis is the. projection 
of a rhombohedral unit cell edge on a plane perpendicular 
to the triad axis. These projections are shown in the 
plan-view diagrams in Figure III j. If the steady applied 
field Hy lies in the x-z plane (¢ = 0°) then the two sites 
are obviously equivalent. However, if ¢ is not equal to zero, 
then the two sodium sites are inequivalent from the point of 
view of quadrupole relaxation. This can also be shown 
theoretically by considering the expressions for Wi and Wo 
in terms of the M-component tensors. These M-components are 
determined by the local environments. Therefore, to have 
identical M-components, in particular M113 and M403" for 
each =P NG Site,the Orientation of the x vaxtseam space for the 
two sites must be different. This would mean therefore that 
for an arbitrary magnetic. field direction the angle would be 
different for the two sodium sites and hence the Wy and W, 
for the 23a nuclei located at these two sites will in general 
be different. 

The single crystal of NaNO, used in our work was in 


the: form of a right circular cylinder, 0.9cm long and 1.0cm 
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Caption Iii. k 


Figure III k is an x-ray photograph taken of our 
Single crystal of NaNO, using a precession X-ray camera. 
The incident X-ray beam was directed along the triad axis 
of the crystal. The crystal was mounted with the cylinder 
axis horizontal. The interpretation of the photograph 


enabling the determination of the orientation parameter $¢ 


is discussed: in the text. 
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in diameter and was obtained from the Harshaw Chemical Company 
Gteveland, Ohio. Ine crystal was cut such that the: triad axie 
was perpendicular to the cylinder axis. The sample was mounted 
in the magnet with the triad axis perpendicular to the vertical 
ax Spel etOrection. .Thus,.byv.,.rotationrgthe.angleswthat, the,triad 
axis makes with the magnetic field could be varied from 


Seton labo. 


0 

As spectroscopic analysis showed the following im- 
PuLitiess: Al-3ppm;.Caz8ppm;, Mg=3ppm;,.Cu<lppm.and.Si<lppm.. .No 
trace of the paramagnetic ions Fe, Mn, and Ni was found. there- 
by indicating that magnetic relaxation should be negligibly 
small at»room. temperature. 

The orientation of the x-axis (the projection of a 
rhombohedral unit cell edge on a plan perpendicular to the 
triad axis) relative to the crystal cylinder axis was found 
using X-rays. Figure III k shows a photograph taken with a 
precession camera with the X-ray beam directed along the triad 
ast For this photograph the crystal cylinder axis was 
horizontaL.(parallel.to Ehe.shegt,edge.of the photograph). 

The diamond pattern of the spots is characteristic of the 
pattern obtained for.a.rhombohedral_unit;cell. s.The,angle 
between the x axis and the cylinder axis of the crystal is 
found tobe 32..8 1.5°. Since the crystal cylinder axis was 
vertical for most of the NMR measurements, the azimuthal angle 
@ was equal to 2.8 1.5°. For later measurements $ was equal 
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* We are grateful to Dr. M. J. Bennett.of the Department of 


Chemistry for taking the X-ray photographs. 
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J) Experimental procedure 


We shall discuss in detail the experimental procedure 


BOtmacncs Particular crystal orlentation @ =12.6 2 71> sand 
6 = 0° where the eylinderical axis of the crystal is vertical. 
Myece “Gases were Considered at this: orientation: (Case I, Ehe 


enhancement of the high frequency satellite due to severe 
Saturation of the centre line. Case II, the enhancement of 
the low frequency satellite due to severe saturation of the 
centre line and Case III,the enhancement of the centre line 
due to severe saturation of the high frequency satellite. 

For Case I measurements, the centre line frequency was 
accurately determined in both the frequency increasing and 
decreasing sweep modes. As described earlier in this section, 
the mean of these two values gives the true centre of the 
centre line. This is the. frequency of the saturating signal, 
the signal generator being kept within + 30Hz of this value 
during- saturation runs. By measuring the satellate frequencies 
at the same time,it was possible to check the correctness of 
the ervystal orientation ©. Also, by monitoring the high tre— 
quency satellite frequency it was possible to detect slight 
frequency changes caused by the slow drift of the magnetic 


field. The saturating) stgqnal could therefore -be adiusted to 
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compensate for this field drift without separately measuring 
the -centre line frequency. The temperature coefficient of 

the quadrupole coupling constant (Andrew et al., 1962) is very 
small at room temperature so one can assume that it is essen- 
tially temperature independent in our laboratory. 

In order to ensure that the centre line is completely 
saturated, we first measured the enhancement of the high fre- 
quency satellites a function of the saturating power BOBS 
applied to the centre line. The observing power applied to 
the high frequency satellite was kept as small as possible 
during these runs to minimize its effect. The experimental 
runs consisted of recording enhanced and unenhanced high fre- 
quency satellite resonances for different values of the satu- 


rating Power P. being applied. A measure of this power was 


at 
provided by the voltage across the saturating coil monitored 
by an oscilloscope connected via a X50 by x50. To improve the 
accuracy, resonances were recorded several times, typically 
three» enhanced and four unenhanced resonances recorded alter- 
nately. Throughout) the enhanced runs, the balance mode cri- 
terion of 'two centimeters' was strictly observed. 

The data is given in Table IV a and is shown in 


Figure IV a, The solid circles refer to measurements made 


with a modulation amplitude of 1.3G (1.46kHz) and the open 
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circles are for measurements made with a modulation ampli- 
tndeoo’: 0265Go(0.¢38kHZ) ait This | datanwassfitted. tor thesBPP 
equation+ (II.65) using a least squares method. The fit is 
shown as a solid line in Figure IV a. The figure clearly 
shows the plateau region that results from the essentially 
complete saturation of the centre line component. From the 
theoretical curve it is found that the enhancement for Baie = 
1600 volts? (corresponding to 40volts peak-to-peak across the 
Saturating coil) differs by less than 0.1% from the 'ideal 


enhancement' with Ey equal to infinity. The extrapolated 


at 
enhancement corresponding to the infinite saturating power is 
Shown on the figure by the dotted line. Accordingly, this 
value of em was selected for all the measurements at this 
crystal orinetation. Furthermore, since an error of 0.1% is 


significant, Ascorrectionsfor- a non-ingtanitenP was always 


sat 
applied. In all cases this was done by measuring the en- 
hancement as a function of Puse and fitting to the BPP equation 
as described. 

For the above results, the spectrometer power Pane 
was-set at a value corresponding to 20 millivolts peak-to-peak 
across the coils. The spectrometer level was indicated by 


the detector level meter. This was calibrated directly in 


terms of the millivolts peak-to-peak across the spectrometer 
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coil by measuring this voltage using an oscilloscope. However, 
the oscilloscope could not be left permanently attached to 
the spectrometer since it acted as an antenna leading to un- 
wanted interference. 

From the data in Figure IV a,we see that the enhance- 
ment is not dependent on the modulation amplitudes of 1.36 
and~0.65G. Asa result: we chose to use a modulation amplitude 


of 0.975G for all the data taken atthe’ crystal’ orientation 


The enhancement of the high frequency satellite was 
now measured (using an amplitude of modulation equal to 
0.975G) as a function of the spectrometer observing power, 


P , with P =~ 16010 funute-. In the case of the lower 


obs at 

observing powers we ran five saturated and seven unsaturated 
resonances. Fewer were required for the higher observing 
powers because of the improved signal-to-noise ratios. Once 
again the balance criterion was strictly adhered to, results 
violating the criterion being rejected. The data is shown in 
Pigure. r¥@b “as "Case “Pt * “The stright *line represents a. feast 
Squares fit or the’ datas” A“statistical”test*was=used*to 
determine whether data at large sabe 3 values should be rejected 


because they were outside of the linear region. The enhance- 


ment values shown in the figure have been corrected because 
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of the non-infinite saturating power as was discussed earlier. 
By substitution of the extrapolated enhancement at zero Pe 
anto equation (11.60). the ratio of the quadnupeleptransztion 
probabilities W/W, may be calculated. 

At this same orientation the above process was re- 
peated for different resonance components. Case II, the en- 
hancement of the low frequency satellite as a result of 
severe saturation of the centre line was measured and also 
Case III, the enhancement of the centre line due to saturation 
of the high satellite was investigated. These two cases are 
Shown sin. Miqune lV. b.wne dala for Case dalvand Case’ Ill as 
given in Tables IV b to f. This data is discussed in Section 
Ly. 

Measurements Similar to Case I were made at crystal 
Orientations Gas. 225. a ef 220 on Tide One t0 A0y Oe Ube, ella eS — 
13550, 15040. 15%.,5 4..and +80. O.sul tisha eee Finally, 
the crystal was remounted in the holder in such a way that $ 
was increased by 30° to 32.8 + 1.5°. This was done bie ro- 
tating the crystal through 90° (this being equivalent to a 
othe change) about an axis parallel to the symmetry axis. 
Measurements corresponding to Case I were then made at this 
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SECTION IV Results and discussion 


We first discuss the enhancement E of a component 
due to severe saturation of another component. The 
enhancement of the population difference is taken to be 
the ratio of the peak-to-peak intensity of the resonances 
with and without saturating power. This is the case 
only if the lineshape is the same in both cases. Such a 
Situation is approximated by the use of small ae values. 
It is essential therefore to extrapolate to this limit 
for quantitative measurements. 


The main purpose was to check the validity of 


applying equation (II 65) to our experiment. The measure- 


O 


ments were carried out in detail at @ = 0 and ¢ = 
2.8 + 1.5°. The data is given in Table IV a and is shown 
in Figure IV a. he error ‘bars are, estimatedgerrors 


determined by the signal-to-noise of the resonances. 

The solid line represents a least squares computer fit 
of the data to (II 65). The good agreement confirms the 
correctness of applying the BPP equation even though 

Pp pealtnoughssmald, iS non-zero. «lie finite ere 


obs 
will have the effect of changing 
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Table IVa 


O O 


= 00 and’ .o-= 248 4 1.5 
The enhancement of tthe high frequency satellite as a function 
of ‘the saturating power applied to theircentre mine. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 


Modulation amplitude is 0.65 G peak-to-peak. 
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ance Enhancement 


(Carb. Un tes) 


14.00 S20 5 O08 
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Modulation amplitude is 1.40 G peak-to-peak 
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Caption: 1Ve‘a 


Figure IV a shows the enhancement of the high 
frequency satellite as a function of the saturating power 
applied to the centre line. The spectrometer observing 
power during these runs corresponded to 20 millivolts peak- 
to-peak measured across the sample coil. The solid line 
represents a least squares fit of equation (II 65) to the 
experimental points. The dotted line is the extrapolated 
enhancement of the satellite for infinite saturating power 


applied to the centre line. 


. Lyi Oe ie 
eft ett 9 4 hemeortaciae ait ore. 6 wa ret 
 ¥awog. oetbteretee: jets aan & er 


a 


-— ace 


efit belos SAT «A EQ0 / tiga’ om sens bowen 4 
ot od (0 (TE) ROS SRB ioe Pres 2% ipe iisenay lll | 
hetsloaatsxs Bye) at ong ‘betseb edt Log tI tn el 

1Swog priterwtse os tntint tot onattatee ods to : | 


co 
2 


~weeg atlovikiim OS Oats fish 


yOl 


SLINN AUNVULIGUV NI #°%d 
zOl l 


(Z7H1EZ0)9690 ‘1dWV ‘GOW d Ol dO 
(ZHI9OVL) OE lL Idwv GOW ‘d Old @ 


LINIWSINVHN&I 


1O-~ 


‘ 


ie 
oe, 
aa 
~~ 
> ’ 
S 
2 ¥ 


; 


10 . 


Oso b WOR’ ¥War’ GO2ELOVIKHT) 


rae@ (3 4erHs)} 


@riob woo vwer 


128 


the ordinate scale however this is of no consequence since we 
are only interested in checking the validity of the BPP equa- 
tion and not in the determining of any absolute enhancement 
values. One can however.use the data to obtain the correction 
factor: for thernen-infinite«saturatingspower, that will be 
applied in the measuring of the enhancement as a function of 

P obs? The dotted line in Figure IV a represents the extra- 
polated enhancement at infinite saturating power. 

Having confirmed the correctness of equation (II 65) 
it was necessary only to take some 7 points on the E/P oat 
curves for the other crystal orientations in order to find the 
saturation correction factor. The E/E oss data for the other 
crystal orientations is presented with the appropriate E/P obs 
data. As a final point we note from Figure IV a that the 
enhancement values are independent of the amplitude modulation 
in contrast with data in Figure IlIe. All the data taken at 
other crystal orientations employed a modulation peak-to-peak 
amplitude of 0.97 G which was adjusted where necessary to 
maintain a fixed modulation amplitude/linewidth ratio. 

The enhancement as a function of Dobe was measured at 
6 = 0° ana @ =h208y% 1.5° for three specific cases: 


Case I, the enhancement of the high frequency satellite due 


to severe saturation of the centre line. Case II, the en- 
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hancement of the low frequency satellite due to the severe 
Saturation of the centre line and Case III, the enhancement 
of the centre line due to the severe saturation of the high 
frequency satellite. 

The data for these three cases is given in Tables IV b,d, and 
f and is shown in Figure IV b. The error bars represent 
estimated errors. It is clear that the enhancement exhibits 
a linear dependence on the small values of the spectrometer 


observing power Po as predicted by Hughes (1966). The 


bs’ 
straight lines in Figure IV b represent a least squares 

computer fit of the data where the weighting factor for each 
point is taken as 1/(estimated error) ~*. By examination of 

the, deviation of each point from the least squares fit, some 
points at large Pops Were rejected as being outside of the 
linear region. 

The extrapolated enhancement at P ech 0 is Significantly 
different.for the,two cases I and, I1.,-The,enhancement.is in- 
deed not expected to be the same in these two cases since the 
energy level spacings are unequal. It was for this reason 
that we developed the rate equations given in Section II taking 
into account this unequal,splitting factor in the hope, that 


our experimental data would be sufficiently accurate to illus- 


trate the effect. Using equations (II 60) and (II 61) we find 
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Table IV b 


6 = 0° and o = 2.8° + 1.5 


0 
Case I, the enhancement of the high frequency satellite due 
to severe saturation of the centre line as a function of the 


spectrometer observing power. 


The conrection»forthe «non=-infinite saturating power is 1.001 
The correction 6 = 0.022 
Pane Unenhanced intensity Enhanced intensity Enhancement 
(aeb.) uni te) (arb. units) (aEb Unies) (corrected) 
225.00 arorh 2) O. 038 6.982 2 0200 Lasoo 2.40. O14 
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2500.00 Tegeses, 0.026 Ha Degg 74 6) 0,032 Lee 2/5 0.007 
3600.00 (EAS PAY 0.020 AB Rasy oy 0025 ie. Ou 0.005 
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10000.00 Byeeileits) O2025 92630 8 pa 9 ay 1.641 0.005 
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Table IV c 


dhe al 


The enhancement of the low frequency satellite as a function 


of the saturating power applied to the centre line. 
spectrometer observing power corresponds to 20 millivolts 
peak-to-peak measured across the sample coil. 


Modulation amplitude is 0.65 G peak-to-peak 


eae 
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Pat Enhancement 


(arb. units) 
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Table IV d 


and ¢ = 2.8 + 1.5° 
Case II, the enhancement of the low frequency satellite due 
to severe saturation of the centre line as a function of the 


spectrometer observing power. 


The correction for the non-infinite saturating power is 1.001 
The correction, 6. =.— 10,022 
ae Unenhanced intensity Enhanced intensity Enhancement 
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Table IV e 


@ = 0° 


and =| 2.0 2 1.5 
The enhancement of-the centre line as a function of the 
saturating power applied to the high frequency sattelite. 

The spectrometer observing power corresponds to 20 millivolts 


peak-to-peak across the sample coil. 


Modulation amplitude is 0.65 G peak-to-peak 
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Table IV £ 


6 = 0° and ¢ = 2.8 + 1.5° 


Case III, the enhancement of the centre line due to severe 
Saturation of the high.fregquency satellite as a function of 


the spectrometer observing power. 


The. correction, for, the. non-infinite saturating power is 1.0017 
the correction. 6 = 0.022 
ae Unenhanced intensity Enhanced intensity Enhancement 
(arb. units) (arbs ,anits) (arb. units) (corrected) 
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Captron iV b 


Figure IV b shows three cases of the enhancement of 
an 2 Fie resonance component in a single crystal of NaNO. due 
to severe saturation of another component as a function of 
the spectrometer observing power eae Case I is the en- 
hancement of the high frequency satellite due to the severe 
saturation of the centre line, Case II is the enhancement 
of the low frequency satellite due to the severe saturation 
of the centre line and Case III is the enhancement of the 
centre line due to the severe saturation of the high frequency 


satellite. Each solid line represents a least squares fit 


of the data to the enhancement equation derived in the text. 
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that W,/W, =s0952 20.015 for Case vr and W/W, = 0.986 
Onuzo for. Case Il. “The error limits for W/W, are obtained 
from the scatter of'the points from the straight. line, and 
are sometimes referred to as external errors. The error 
limits on W,/W, as obtained from the estimated errors are 
invariably larger than the external errors usually by a 
BACzLOrsCfad OF S.aynthissis due*to the difficulty of 
initially estimating the error limits in the measuring of 
the resonance intensities. For Case III, the value of 
W,/W, derived by extrapolation is 0.988 + 0.022, in 
agreement with the values obtained for Case I and II. 

Such agreement for W,/W, would not be expected. if 
Magnetic relaxation due to paramagnetic ions were signif- 
icant (BPP, 1948). We denote magnetic transition proba- 


bility by (Andrew and Swanson, 1960) 


W =, ex 0i+om) / (Lec mo tel) W, CEVe 51.) 
m>m- 1 
and assume W, to be same throughout the. crystal and the 
same for each transition. (This will be true for direct 


relaxation. by a paramagnetic ion but net +£ the relaxation 
were determined by a spin diffusion). The rate equations 


for Case I, where y = W/W), are 


N/Wi= BE le xtoy) aN, (2xr8y tere aN hla Xt Sy) gan eke a! 
(IV 2) 
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Ny = “Ny (2+x+6y)+4N (yt+P,)+N_) (x) +2n, [l+y+6 (1+3y+x) ] (IV) 3) 
Ny + Ny (x) +4N) (y+P.)-N_, (xt2+6y)+2n. [l+y-6 (1+3y+x) ]. (EVI4) 


Solving these equations in the usual way we find that the 
Case I enhancement, for zero observing power, is given to 


PieSt order im yiand © by 


E = [1+2x] [1- 3xy -@2 fxdoi (IV 5) 
1+x CEE Ce a ee Ree | 


The enhancement for Case II is from symmetry arguments given 


by 


E = [1+2x 1- 3xy + xd | (IV 6) 
1+x (1+x) (1+2x) 1+2x | 


For Case III the rate equations are given by 


No = N, (1-x+3y+3P,) “N, (2x+8y) +N_, (1-x+3y)+2n) (2xty-1) (2V87) 
ae aes 

Ny N, (2+x+6y+6P,)+4N y+N_ x + 2n,{ity+s (1 3yt+x) ] (IV 8) 

No. = = 3 ry-6 (1+3 IV 9) 
Ny N,x+4N oy N_, (xt2t6y) + én, [ity ( y+x)], ( 


and the enhancement of the centre line for zero observing 


power 26 given to first sorder yim y and =o by 
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(TV 10) 


E = }|2x+3 i+ 6 (14x) - y (7x7 4+8x44) 
x+2 2x+3 > lp aars) OES PCy i 


Substituting the measured enhancement values for the Case I, 
Pend tm srmto ‘equations, (my 5), (IV 6), andi (ay Oe we 


find using a least squares analysis that 


y = W,/W, = -0.00404 0.0034 (IV 11) 


This shows that magnetic relaxation effects are completely 
negligible in our crystal. This is to be expected from the 
spectroscopic analysis given in Section III 8 where no trace 
of the paramagnetic impurities is indicated. 

The enhancement data obtained for the various values 
of 89 at ¢ = 2.8° are given .in Tables, .EVije +0299 Further 
enhancement data obtained when the crystal was tipped through 
90° as described in Section III are given for 9 = 32.8° in 
Tablessa, 6, yauCpOe, 4, G, and c. <I8sem thévemhancement data, 
the values of W/W, were obtained as described above, and are 
lzsted@in Table iV 'T fon -o)= 2.8° and in Table IV = for 
oc= 3986. The data for od = Dee. is plotted.in Figuresiv ¢ 
anduter consistency the value at 6 = 0° is that obtained for 
Case f, since Case Il and Case Iilewere only carried out at 
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po= 0.2) Theserror lamits. The data for o7 S228" sis plotted 


in Figure IV qd. The error limits shown in these figures are 
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Pape. FV «¢ 


@ = 22.5° and ¢ = 2.8 + 1.5 
The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peka-to-peak measured across the sample coil. 
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Table IV h 


3 22.508 e ¢ 2-2.¢9 4-2.5° 
Gase “I; “the “enhancement of the high frequency satellite’ due 
to severe saturation of the centre line as a function of the 


spectrometer observing power. 


The’*correction for the non-infinite saturating power is 1.0013 
The correction 6 = 0.017 
Ee Unenhanced intensity Enhanced intensity Enhancement 
(arb... units) (arb. units) (arb. units) (corrected) 
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Table IV i 


@ = 42° and 6 = 2.8 +.1.5 
The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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Table IV j 


O O 


and @-2 2,8 2 1.5 
Case I, the enhancement of the high frequency satellite due 
to severe saturation of the centre line as a function of the 
spectrometer observing power. 


The correction for the non-infinite saturating power is 1.001 


The correction 6 = 0.007 
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Table IV k 


6 = 74° and o = 2.8 + 1.5° 


The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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(arb. Units) 


26,010. 1:20 5) O05 
5308 fl Sao Se hf, 20S 
fle) 1.454 O2OS 
FO. O0@ 1.464 ONO'S 
20°00 era S Ores 


40.00 EE473 OS 


~~ 
~- ae 
< 
p a 
erescsy F 
BS i“ 
* >> 


qeerres 1atog | pritvasado -eiteind $9 eye 


ott sxotos “Bowwesom “ars 3-48 


i: 
| ee 
7 
a, ee 
. 
i ee 
1 
| . 
iB ibaa 
& oe 
¥ 74 


Op ly 4 %. ie = ey ‘bis cs TS 


aye 
fai 


a 


A sige oka As Yoomebanit 


ate 


| ie afi 
bontqgs tewog paitenbdse” Sif 


aoa it 


, 1 ‘ ; 7. 
= 7 
ee: 4 


(ation .dxs) 


é 
ee a ia aon Ae 
meee) ae 


HEALY tithin 


145 


Table IV 1 


6 = 74° and » = 2.8 + 1.5° 
Case I, the enhancement of the high frequency satellite due 
to severe Saturation of the centre line as a function of the 


spectrometer observing power. 


the, correction,.form.the .non-infinite, saturating power is 1.001 
The correction 6 = 0.008 
Aisa Unenhanced intensity Enhanced intensity Enhancement 
(arb. Units) (arb. units) (arb. units) (corrected) 
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Table IV m 


6 = 90° 


anasO., 22652 195 
The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 


spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil 
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Table IV n 


6 = 90° and # = 2.8 + 1.5 
Case I, the enhancement of the high frequency satellite due 
to severe saturation of the centre line as a function of the 


Spectrometer observing power. 


The correction. for the.non-infinite saturating power is 1.0007 
Ene-correction ¢6 ="0.01L1 
Pane Unenhanced intensity Enhanced intensity Enhancement 
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Table IV 0 


f= 100° and ¢ = 2.8) #°1.5° 

The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak across the sample coil. 
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Table IV p 


@ = 100;sand ¢.= 2.8 + 1.5° 

Case I, the enhancement of the high frequency satellite due 
LOpecveres.sakurationuot the,centre,line.as a,function ofthe 
spectrometer observing power. 


TieRecOLrect ion itor. thesnon-infinite saturating) power 1s 1.0015 


DHescorrection: 6 = 0.010 
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Tabie IV q 


© = 112.5° and ¢ = 2.8 + 1.5 
The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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Table IV r 


@ = 112.5° and # = 2.8 + 1.5 
Case I, the enhancement of the high frequency satellite due 
toe Severe .saturation,o£f’ the, centre, line as\iavfunction “of: the 
spectrometer observing power. 


Thencorrection, fon the non-infdind teysatunating power: is. 1.0009 


The correction 6 = 0.006 

Beis Unenhanced intensity Enhanced intensity Enhancement 
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Table IV 5s 


6 = 135° and o 2.8 + 1.5° 


The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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Tabike YEV «t 


@ = 135° and » = 2.8 + 1.5° 
Case I, the enhancement of the high frequency satellite due 
to iseveresaturation;,of. the: centre slbine.as a functionmiof the 


spectrometer observing power. 


the. -correctiom fon the, non-infinite,saturating power is 1.0013 
The correction 6 = 0.006 
Pee Unenhanced intensity Enhanced intensity Enhancement 
(arb. Units) (arb. units) (arb. units) (corrected) 
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Table IV u 


6 = 150° and @ = 2.8 + 1.5° 

The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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Table IV v 


& = 150° sand o.= 2.8 + 225° 


Case I, the enhaneement of the high frequency satellite due 
tm severe saturation. of the. icentnre. line .as.-a--function sof, the 


spectrometer observing power. 


The correction for the non-infinite saturating power is 1.0009 
The correction 6 = 0.014 
P aie Unenhanced intensity Enhanced intensity Enhancement 
(ane. Uli S) Carb units) Carb’. units) (corrected) 
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Table IV w 


B= ts7°o 


and Ora Oa Le 5 
The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 


spectrometer observing power corresponds to 20 millivolts 


peak-to-peak across the sample coil. 
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Table IV x 


@ = 157.5° and ¢ = 2.8 + 1.5° 


Case I, the enhancement of the high frequency satellite due 
to severe saturation of the centre line as a function of the 


spectrometer observing power. 


The correction for the non-infinite saturating power is 1.0016 
The correction 6. = 0.017 
P obs Unenhanced intensity Enhanced intensity Enhancement 
(arb. units) (aro. Units) (arb. UnEts) (corrected) 
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Table IV y 


pe =) 1800 ands o.=02.8 + 4.5° 


The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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(arb. units) 
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Table IV z 


6 = 180° and ¢ = 2.8 + 1.5 
Case I, the enhancement of the high frequency satellite due 
to severe saturation of the centre line as a function of the 


spectrometer observing power. 


thes correction for’ the*nen-infinite’ saturating power is 1.0007 
Pae correction 6 = 0.022 
Beas Unenhanced intensity Enhanced intensity Enhancement 
(arb. units) (arb. Units) (arb. units) (corrected) 
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Table IV = 


@ = 42° and $ = 32.8 + 1.5° 

The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
Spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 
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sat 
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Table IV 8B 


of 4209 was = 328.3 4.52 


Case I, the enhancement of the high frequency satellite due 
toethe severe saturation of ‘the: centre! line’ 'as® a. function® of 
the spectrometer observing power. 


The correction for the non-infinite saturating power is 1.0010 


The correction 6 = 0.008 
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Table IV y 


6 = 69° and ¢ = 32.8 + 1.5 

The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak across the sample coil. 
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Table IV 6 


O O 


and Or= 3u56 = 1.5 
Case I, the enhancement of the high frequency satellite due 
to the severe saturation of the centre line as a function of 
the spectrometer observing power. 


The correction for the non-infinite saturating power is 1.0016 


The correction 6 = 0.007 
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Table IV ¢€ 


Gn=) 100°. ands du=23268-49155° 


The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer observing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 


P Enhancement 
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Table IV n 


6 = wl 104 Candid s= 43208 et 01 25° 


Case I, the enhancement of the high frequency satellite due 
tothe wsevere saturation of, the.centre line .as a function of 
the spectrometer observing power. 


The correction for the non-infinite saturating power is 1.0007 


Maevcerrection-6 =-0.007 

eS Unenhanced intensity Enhanced intensity Enhancement 
(arbe unats) (arb nits’) (arbS -uni-es) (corrected) 
400.00 2.867": 02038 4 S2h382 0-038 LAGS = 10) 5074 
625.00 34250 021050 4-778 0.058 L470 0.029 
900.00 3 £500 0.030 5.179 0.035 L74e0 0.016 
P22 520.0 38675 0203.0 5.454 0.035 eee 020.05 
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2500200 32475 0026 5229 2032 1753505 O05 
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Table IV & 


fem oc Mene’o. = 32:3 = 15° 


The enhancement of the high frequency satellite as a function 
of the saturating power applied to the centre line. The 
spectrometer obServing power corresponds to 20 millivolts 


peak-to-peak measured across the sample coil. 


P Enhancement 
sat 


(arb. Units) 


300 L300 £00003 
5.00 W388 01.03 
7.00 1.463 0.03 
L000 1.463 04,03 
20.00 1.463 OaGs 
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Table IV Z 


6 = 038° tandys =932.8 #1.5° 


Case I, the enhancement of the high frequency satellite due 
to the severe saturation of the centre line as a function of 


the spectrometer observing power. 


Hnencorrection for the non-infinite saturating power is 1.0011 
The correction 6 = 0.007 
Pec Unenhanced intensity Enhanced intensity Enhancement 
(arb. unites) (arb. units) (arb. units) (corrected) 
400.00 Dei 200138 3, Ode O08 1, 490) 240) 032 
900.00 Zi Sel 0 0.030 3 Jao 0035 1.506 02023 
225.00 2 to 0.030 Ae 155 02 035 Lea: Q2021 
1600.00 2 Bw) 2025 4.180 0. 0.29 1.520 G07 
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3600.00 2330 0.020 BTR, 0.025 L5G" 0.017 
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Table IV T 


The data being measured under the conditions of Case I, the 
enhancement behaviour of the high frequency satellite due to 


severe saturation of the centre line 


6 Noa 
(degrees) 

0.0 OF 952. O27 005 
PLUS ORS OO17 
42.0 1.026 0.016 
74.0 0. 29¢ 0.015 
90.0 0). 8:35 0022 

100.0 0.994 O07 
PARTS Le. LL Or025 
13550 Le LoS ero) aky) 
150)./0 Lg 7 OF0.568 
TES Reyes 0.956 07022 
180.0 0.936 0. O16 


6 = 0° for Case II measurements gives W/W, = 0.986 + 0.029 
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6 = 0 for Case III measurements gives Wo/W, = 90938 20.22 
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Caption IV c 


Figure IV c shows the orientation dependence of our 
values of W,/Wy Mi 0 Sea 2 6 Oat 1.5°. The solid line is a 


least squares fit of the Pietila theoretical equations. 
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Table IV 


$= 32. 0%yaeus” 
The data being measured under the conditions of Case I, the 
enhancement of the high frequency satellite as a result of 


the. severe saturation of the centre line 


6 No 
(degrees) 
42.0 1.249) & 02023 
69.0 1.160 0.027 
ELO%.0 0.7871 07-010 
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Caption Gavi 


Figure IV d shows the orientation dependence of our 
values of W,/W, at ¢ = 32.8 #1.5°, The solid line is a 
least squares fit of the Pietila theoretical equations. 


The dashed line shows the least squares fit for > = 2.8 + 


fle ae 


a 


tue) Xo ‘eanetieg PHY 1G 
BBL, anit ‘eine 


0 30 60 90 120 150 180 
ANGLE @ (DEGREES) 


re 


the externalwerrors., obtained from the linear fit°of the ene 
hancement data. This means that our estimated errors are not 
mirectly anvolved in the limits except in so!tar sac they 
determine the weighting factor applied to éach point. Since 
various checks mentioned in Section III showed no evidence 

of non-linearity and since the enhancement values are the 
Simple ratio of successive resonances we are confident that 
systematic errors are negligible in our measurements of 
Wo/W,- We therefore place a high degree of reliance on the 
error limits on the values of W/W, shown in Figure IV c and 


° has been published by Hughes, 


Pigure, LVid.- Theewoerk at_@-=,0 
Reed and Snyder (1970). 

We now consider the comparison of this data with the 
theoretical orientation dependence as developed in Section 
II. Since we have not measured Wa and Wo independently, we 
are unable to determine all the parameters Miia’ Mj 3137 


and M We therefore normalize with respect 


M3333) “Wades hoe 
to M and define the following quantities 
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M23 7 ago ii1: (iy 12 


Assuming that the values of the M components are the same for 


Wa and Wor the orientation dependence is therefore, from 


equation (II 86) 


e 2 
(1+4M) 333 + 2M3333) + (6- 6M33433) Cos 6 
a = +(1-4M- +2M- Neos ofa he cos3o¢+M; Bn oir anes 
iw ioee 3333 (bea) ele? 3 
2 
i 2 
(444M) 3137M3333) ~ (12M) ,,3-9M3333)c08 8 
i ue 
-(4- 16M) 343+ 8M 3333) -16 (Mj 4 1 3COS364+M, , 538in3>) sin 8cosé 


lek 3) 


The solid line-in Figure. IV c represents.a least square fit 
ef the crystal data to equation (IV 13) giving the following 


Mo.values 
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The agreement between the experimental values and the 
theory is very good, thus confirming the correctness of equa- 
tion (IV 13) and the assumption that the M-components are the 


same for Wi and Wo. As was discuSsed in Section II this 


implies that the quadrupole relaxation is via an indirect 
mechanism. This is not unexpected, since the direct process 
should not contribute significantly compared to the indirect 
process. Moreover we would not expect rapid reorientation of 
the NO groups at a temperature so far below the melting point 


Cf the.crystal. 
Figure Iv d shows the W/W, data fon Ovee32700= RS oO. 


The purpose of these measurements was to determine M: cosh + 


LES 
M3 4238in¢ from a value of different from that in the first 
case so that by combining with the data at > = 2 >a one can 


find M123 and M3413° 


fo) 
given earlier, we found at @ = 32.8 :1.5°, 


Using the values of M3137 M3333 


cos3¢ + M: sinsé = 07162 = .0.014 Ev i) 
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Prom {iy 16), (LV 17) we ftind 


S = -0. 20 20ks 
Mi 113 Obs 


and GLY 713) 


% = =O) 
M1123 0.147 
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These runs were made at 0 = 42°, 69°, 110° and 138 where the 


term sin’0cos® is relatively large. This enables M3 413C°8¢ + 
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Table IV A 


M) 313 = 0.85 + 0.015 

M3333 = /0.822 + 0.016 

M1 4730083 + Mj 19382139 = -0.092 + 0.014 
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My1238in¢ to be determined with the most precision. Values 
between 42° and 69° and between 110° and 138° were avoided due 
to the proximity of the satellites and centre line. 

We now discuss results of the orientation dependence 
of Ws and Wy published by Niemela (1967) during the course of 


this work. Niemela measured Wo and Wy for Are in a single 


crystal of NaNO, at 77°K for several values of 6 in the range 


9° 


to 90°. The data is shown in Figure IV e. In these 
measurements the time dependent behaviour of the spin system 
following the saturation of a different transition was studied 
using a pulse method. This method was first suggested by 
Goldburg-W%1959) but ds,different from: the wnes achuallysused?by 
Goldburg. Niemela measured the resonance intensities as a 
function of the time after applying the saturating signal. In 
contrast, Goldburg only measured the resonance intensity imme- 
iately following the application of the saturating signal. 

As was mentioned in Section I, we did not use a pulse 
method because of the danger of exciting other components 
with the 90°pulse. Strictly speaking,the duration of the 90° 
pulse should be much longer than To, the spin-spin relaxation 
time. Moreover, cee should be much less than the separation 
between the adjacent 2 ie resonances. The pulse duration 


’ -4 
used by Niemela was 40 microseconds. Since T,, is~l0 seconds 


er ae a tN ve cla 


~ 
ss 
“ve 
i 
ion 


BYVi 


eae V motereotg Jeon ad itd | / oc OF * 
seb Bebiove exow "SEL baa. "OEE ‘nsowded ons te bas | sh a 
.ofil sxtiso. bas est ilietsa ‘of? 20 yeulsony ia 

oneboedeb notietteito ems to ativast sev0dib We: an 

io Betveo edt parseb TARE) alomeid ya bedettdua re bas 
= trike 5 oe apt xto% iv bn eee: hue << 
apres ons Gi § to sattey sien iets tof x We: ds te tot ‘ 
gaent mi .a VI one BE “twodte: ‘ak ageb edt S00 of 

moteya mige ent to mmokvsded inpbnogeb, emit esis 
batbuse esw sotstedar tosx@R2ib 6 to aotter388 ods 
ve hetaspous sexta Baw Bodsom aif? Boston setug 6 pe 

vd bees yileugos| sho edd make gnexe232b ei sue (eel) px 


S as soltiemstnt somanoees orld boxwasem siomoin 
at .fsapbe patdexrwdee odd pniylage x09%e oma oon to 


“om. yd ad et Sonenoeee, add boxesem bLoo 
copie pniderdne elit #0 wotteok i Ne enone 

selva s Sey tom bub ow vt woldpes fal an tas: @ es GA | . 
einenogmon x16lfo paid tone to nepash os 0 eassosd bortsam 

“02 Sid Yo nottewyb. orld epritieods yisobrse “ealug?oe. on iidiw 
foitsxsiox niqa-aige aug | vg? neds teprol. ea ed. Biuede ealug 
ntokisiedse oft asd est doum ed biuode * stevoex0M vom 
Hortsieb tes ont .sapnenoRet ee snsostbe sit nsewsed: ay 
abmoose O imai gt eomte -ebaogesex9.im Ob a sfomoi xd bean: 


° 


¥ 


a 


bis & a 
ae p aa 
. 


: : Y « ; i 
| aS i p 
7 _ 


177 


Candy rs soa) 


and the minimum resonance separation (at 6 = 40 
was approximately 60 kHz away from the central value (8MHz) , 
it is clear therefore that Niemela was in fact observing an 
induction decay signal which included a significant contri- 
bution from a component other than that being studied. Since 
the recovery behaviour of all the resonances is not the same, 
this admixture of resonances will give rise to errors in the 
measured values of Wi and Wo, 


It should be pointed out that the magnetic relaxation contri- 


especially at 0 = 40° and 73.5°. 


bution in Niemelas crystal at 77°K was by no means negligible, 
the value of W/W, being approximately 0.10, at least 30 times 
larger than in our crystal. While it is possible to correct 
for the effect of magnetic relaxation it is necessary to assume 


that, the W, is unztform throughout the sample. For large 


3 
magnetic effects, this may not be valid. 

In fitting the data to the orientation dependence 
given by Pietila (1968), Niemela_used an incorrect symmetry 
argument in showing that M) 11300839 + My 19381n3¢ LS "Zero. 
According to Niemela the symmetry of the NaNo, crystal 2s. suach 
that the crystal is invarient under the rotation 679 i+ ™ and 
therefore M, 1130083? + My 42351n3¢ as given in equation (II 86) 


must vanish. This is incorrect for two reasons. Firstly, 


the NaNO. crystal is not invariant under the rotation 078 + T. 
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Secondly, even if it were,it does not follow from equation 
(20.86) that the term M)413C083¢ + Mj 49381n3¢ is zero. The 
dashed line in Figure IV e represent a least squares fit of 


Niemela's data to equation (II 86) with M sin3¢ + M 


TES 1128 
cos3$)| set" to’ zere, “This fit differs slightly from'Niemela's 


because of a partial error in his theoretical expressions 


(Pietila, 1968). The corresponding M-values are given by 
re 4 -2 -1 
e'Q Miill ea Se 0 LO seconds (IV -19) 
22 a -2 =a 
e'Q M313 = seo 009% a 0 seconds (TV 20) 
Be2 =2 -] 
e'Q M3333 = 512029 x0 seconds CTW 2 15) 


where the error limits are external errors determined by the 
scatter of the data points about the theoretical curves. 
While the fit between experiment and theory agree within the 
error limits given- by Niemela, the fit is unsatisfactory 
because most of the Wy points lie below the theoretical 
curve and most of the Wo points. lie above the theoretical 
curve. We have therefore re-analyzed the data with M413 


cos3¢ + M sin3¢ # 0. A least squares fit of the data is 


1123 
ehown as a solid line in Figure IV e. For the crystal ori- 
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entation used by Niemela, the angle ¢ is 0° when measured with 
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Caption IV e 


Figure IV e shows the orientation dependence of Wo 


and W. obtained by Niemela at 77°K. The dashed line 


al 
represents a least squares fit of the data to the theoretical 
form given by Pietila with M413 = 0.. The solid line 


represents a least squares fit of the data with M413 ~ 0. 
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respect to the projection of a rhombohedral unit cell edge in 
a plane perpendicular to the symmetry axis. In such a co- 


ordinate system, the M-values are given by 


6-0°M, 444 = 2.42 + °0.05 X100- seconds — (FV 22) 
e7Q°M, 413 = 1-84 + 0.05 x10” seconds * (fv 23) 
SO aMea ad =68.& 10.05) S10 Isecondame (IV 24) 
BGO Mant. POS tae Wo. x0 eecondes se (7.95) 


The importance of the M1313 term is more clearly 
illustrated in Figure IV f where the ratio W/W, is shown as 
a function of 80. The best fit with Mji13 7 0 is shown as a 


dashed line and the corresponding ratios of the M-values are 


Meaa/) yn ae? eyo bo (IV 26) 
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Caption. IV £ 


Figure IV £ shows the orientation dependence of the 
ratio W/W, obtained by Niemela at 77°K. The dashed line is 
a least squares fit of the theoretical form with M4113 = 0. 
The solid line is a least squares fit with M1413 720. 
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M3 333/ Miqq1 = 0.66 + 0.04 (IV 29) 


= 0.26 


lie 
oO 
e 
(=) 
NO 
e 


Pie, yi (IV 30) 


We conclude that the experimental data are much more 
precise than was first believed (Niemela, 1967), the a- 
greement being highly satisfactory so long as M) 1130083¢ + 
My 4938in3¢ is not set to zero. 

We now compare the experimental values of the M- 
components obtained at TT ER by Niemela (as re-analyzed by us, 
Hughes, Reed and Snyder, 1970) and by ourselves at room tem- 
perature with the theoretical predictions given in Section 

O We 
£P, SOR. 727K, Mijq1° M) 313 and M3333 are all positive as ex 
pected for components. of the form Miu gaB° Our values are in 
agreement with this since Mi 313/ Miii1 and M3333/ Miia are 

a 2 (ope 
also positive. Also, M3413 at 77 K is less than My 412"1313° 


, 2 
Similarly at room temperature we find that My3135 M1111"1313 


2 
and M1123 < My 311!2323 (since it was shown by Snyder and Hughes 


(L970). that M5393 = M,313)° 


We now consider the relationship between wy and Wey 
the transition probabilities for the spin systems A and B and 
also the relationship between WS and Wo For simplicity we 


first suppose that the nearest neighbour oxygens alone are 
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responsible for the time-dependent electric field gradient 

ae the particular - ne Site. This» assumption is not unrea-— 
sonable since as will be seen the quadrupole spin-lattice 

8 =a) 

Yr : 


interaction falls off very rapidly as: r ~ or The 


relationship of these six nearest-neighbour oxygen atoms to 
the ee site is Shown.in Figure il 4 forvcases A anes. 

As was pointed out in Section III 8, the two lattice sites are 
equivalent apart from a translational and rotational oper- 


fe 


ation. It is convenient to define a coordinate system (x”’, 


4- Figs 2 


Vo. GZ for the spin systems A lattice site where z°~ axis 


coincides with the triad axis and the x~~ axis is orientated 
so thet two oxygen atoms lie in the ™ =z “-plane, “ftiis clear 
from Figure IIIj that this x°°z°~ plane (there are of course 
3 of them orientated at 120° to each other) is a mirror plane 
for these nearest neighbours (however this is not so if one 
considers further neighbours). We denote the angle between 
the x~* axis and the x axis (defined previously as the pro- 
j7@ction of the rhombohedral unit cell edge on a, plane per- 
pendicular to the crystal triad axis) by do: 

A plan view of the. situation is shown, in Figure.IV g. 
From the crystallographic information the angle o,is calcu- 
lated to equal 18°. For the B lattice sites one can define 


CA 


° AAs Are 2Ae 4 
a Similar system of axis (x guny. poe such “that. x 


Est 


tnatbere bist? otxaoele ineba 
“BSta ton BL, motsqauees ebaT ste wif sages 
eb izytel-ni¢e slogetSsup ont asea od tiw as conte al 
wet Py ao BE ee yibéqes yxawi to: atts? nolsoa 
od gmOgs agpEKo LU OTHE ERAT OSMAOTE ule eae: e's 
2) bap A 29865 so? -¢ PEL, onupet oa avae at, 0288 
ets setie odittsl ows sdt .8 {11 Apiiteee me, 390, besniog « 
-yeqo lanoltetos bas (shotsalensay 6 moxdadisqs cae £9} 
, x) Metevs ots ibtoen 6 smiteb oF JaeLnevaso at FF ; of ob: : 
eiss ~~ s sisdw stie costaet & ameteye jaige odd wok, Cos | 
botstneix0 si aixs ~ x eae, bas BkxS Rakes eds dtiw, s 
tselo 21 43 a ima: Da oe ods at wit, emoze moons salt 
eamios To Sis exsit) onghe ~ > lad aids sede orer ie 
ensigq tontim 5s et (lise aoe oF Poss, te-bedsdneixo 1 
ado tice ton ef eine xevewor) “eoodeton. Jasssor ¢ 
hsewted of eis ait stone ow  (axpoddpisa redsind . | 
| -O4qg @dd fe wlesoiveng peatteb) bead x fit bas ains xs = +: ! 


- 
=16q snelg s mo Sens ifs stan ere add. 20 ical 


+e VI ae bial Py a ai ons 56 <n nala & ayy ay mtd 
~wolso. ah .¢ aces out nolsamiotas poaaraehica oi fe” 
ae em pe shoe RE ae 3. «SL Loupe os Pedak | 


“ws get dowe (og (Po hy (ot res to metaye islimia 5 


A 


184 


Caption IV-g 


Figure IV g shows schematically the two spin sites 


with their appropriate principal axis and the angles d>5° 
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eae 


Zz plane is again a mirror plane containing two nearest 


eaAesA 


neighbour oxygen atoms. Again the angle between the x and 


a 


van O e 
the x axis is 18°, the sense of rotation from x to x being 


Opposite to that from x to. x ~~ 

The two lattice sites A and B are completely equivalent 
if A and B are viewed from their respective coordinate systems. 
The values of the M-components will therefore be the same for 
the two systems if they are defined in the appropriate co- 
ordinate system. We denote the M-components at A sites that 
are non-zero for the symmetry by Pe M3313" M3333) M3313" 


and M> We can therefore write quite generally 


1123° 


eae pag 3333 


a 2 se, 4eA —mM77 a Z 
Ww = e“Q a +4M gga Sy =7i(12Ms A OMe) eoare 


-(4M3331 16M; 313+ +8M M3333) COS 8 16 (M} 3 13¢083¢ = Mji03 


sin3o~~) Sac e cosé | (IV 3) 


= J 4 A 
Where $°* is the azimuthal angle of H, relative to the (x”’, 


ze 
2- 


y°7, 2°7) system. For our choice of x° axis the x — z 
plane.is a mirror plane if we consider only the nearest. 
neighbour oxygens. The transition probability Wy will there- 
fore be invarient with respect to mirror reflections in the 


x°+z°* plane. Such a mirror reflection corresponds to the 
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Ghange: 0° “+,=0" **in,equation- {IV 3). This invariance property 
requires that M3123 be equal to zero. This can also be shown 
from the invariance of We. 


We denote the M-components for the type B system 


4A 2A 2A 


relative to the (x7 y ) coordinate system by M 


Pe a) 


UR Ei 


Meo Me. 5 M4 --swand Mags Using argument similar to that 


Be 34. 3¢ 3388"%55 1413 Lb235 
applied to type A it follows that M; 


A 


1123 


Since the sites A and B are equivalent when viewed from their 


is zero. Moreover, 


appropriate (principal) coordinates systems, the values of the 


Messe 


Cerresponding M-components are equal. Putting Mii. = M3111 


etc., the transition probabilities can be written 
We aor (AMC AM ees oo) (LM -9MZ2,.,) cos“6 
i. wise tag, ene ea ee 1313 ga 38 


=(4M" -. =16M- 5: .+8M cos *6- 16M 


ibaa ball enhen 3333 


ae ° Ss 
1113¢°83¢ Sin’ @cos@ | 


(IV 32) 


We = (gue [ons jaatMiSa3t? M3333) + (6M) 111~6M3333) C08 0 
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Lig ie 13034 M3333 


yee 
1113°°83¢ Sin @ cosé@ 


(IV433) 


“NOE eee a A ee 


eh ie a ‘- =4 ing is ra tees ap 
yoteqoxg | soaniaswa aia? (2 VI), mo i? | 
Pipes ‘wamat vito we Sen Po 
nworfe ad ents 80 eld O98 ce aes sc 


j ee: ae. Pe dh, 
FT| iy ie Fr we dae ie 


« ~ “As aa A 1 7 
ap ya mosaye at anithroas ( | 


dads, oF xsLimte Jasmupts paied ? ; pe bas - rr 
“SE ML 


| tevosxoM ORgR ex tenth shen 9 2wWO 


shonS megs bowsty -niedw snglaviupe pe * 6 bas~A mete 21 Ss 


| és 20. 200! hey ont vemotuve pedanibsccs ‘Usqioaksa) Ox: 
Gm = rena eoees ae ot6 peck f: 10a: 31 
if ab? 
y 7 or 
istatiuw od mo aoi3 LL idndoxg ota ienexd an a 


5> me 
OP eoSle pee ‘saa 


+. =4ms) es 
i can pee : 
| r o> £ cotta. |e, ik i roe Ry | a if | 428s ; > fh 

| [ Saas ate ge e095 f pp ene 200 ie Pree 

4S vi) ae. ck a 


fs g. A ae ~s TP ‘ 
anes aul 8 ta ha 


bs ote shies wt re 


p4a9 nae Ah. % Ee 


ai 


4 
; 


fv Ls 
ATs ms 
, A ve ny, 
7 waa 2 rag 
fi pinay Da mae ae 
> : CP 4 


187 


W = * 22-2 2ae 
ik as jessy FAM 7 3137M@3333) ~ (12M) 31379M3533) cos’ 0 
mM Phe 2) eM = 21) Gos b= WeMo1 4 cons cass neo eose 
ital 131379!3333 ke peers 
(Iv 34) 
Be 2 2 4re 4Are Are rres 2 
we a [ottreccerey 2M3333) + (6M)}},—-6M3333)cos' 6 


+(M> <7 AM 


4 7 a ZAP ° 3 
11117 1!1313+2M3333) C08 8 + 4M cos3¢°“ “sin” 6cosé | 
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(Gs BY) 


The quantities Wy and Wo, whose ratio is measured in 
our double resonance experiment, are equal to (we + W,)/2 
and (we a W5)/2. Therefore, from equations “(1V.39)) (iv 33), 


Cav 34) ¥and )1V735)""*it: follows ythat 
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But from Figure IV g it is clear that 

Py Sikes o. (IV 38) 
and 

Gy ig OM ao. OMe EY) 


Substituting in equations (IV 36) and (IV 37) we get 


W. = Eo. 


(12M**, ,-9M 
a 48 


ee uw? ? ce 2a 2 
(Corre +4M) 3737M3333) 131373333) Cos 8 


ee 2A se, AeA ° 3 
~ (4M5 34 716M) 31 3+8M3333) - 16M) 113 (cos3¢cos3$) sin 8cosé | 
CEV 420) 


re 4-4 2 
| aaa +4M3 313%2M3333) + (6M))11—-6M3333) Cos 8 
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(IV 41) 


The expressions for Wy and W, are comprised of independent 
M-components only one of which, M3313" involves the angle 
¢°*. Since a redefining of the x°° axis is equivalent to a 
change in the angle ¢°° it follows that the M-components 


and M; will be independent of the choice of 
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the x°° axis. A change in $°* will however affect the M3313 


term. The M3313 term will change in such a manner as to keep 
the value of Wy and Wo invariant fOr the Change in 0 =. "We 


Can therefore write 
Piry Ga: | eee 
JOISTS Bg bees (IV 42) 


SEC ia eee 


where the unprimed M-components refer to the unprimed co- 
ordinate system where the x axis is along the projection of 


the rhombohadral cell edge in a plane perpendicular to the 


4, 


triad axis. Finally it is convenient to replace M}713C083%, 
by Mi34413° Equations (IV 40) and (IV 41) can now be written 
ea (4M +4M -M ) = (2M -9M eee 
Vera 8 jae” gajaen 3683 isdopens ees 
o 3 
ES = 8 iV 43 
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showing that the spin system comprising type A and type B are 
equivalent to a single spin species with M4137 M371 13C°83¢,- 
In the above calculation it was assumed that the relaxation 
was caused by nearest neighbour oxygen atoms. However, the 
result is independent of this assumption as can be seen from 
the following argument. Suppose that the nearest neighbour 
oxygens do not contribute at all to the quadrupole relaxation. 
Then the A and B spin species are no longer different, the xz 
plane is now a mirror plane and ¢ = $¢° “= $°°*. We would then 
get the same results as above except that oe is 0° instead 

of 18°. In reality, the quadrupole relaxation field will be 
partly due to the oxygen and nitrogen atoms and partly due to 
the rest of the lattice neighbours. In considering we note 


Ze se 


that the principal coordinate systems (x ““y~"z"~) and (x°°° 
y°°°z~~*) can always be selected such that M3193 and M3353 are 
zero. Furthermore, from the symmetry of the NaNO. crystal 


2e 


structure, the x~* and x°~~ will be symmetrically disposed 
with respect to the x axis. While it is not obvious what the 
angle between the x and x~* or x and x°~~ axes will be,it 

seems likely that uA should lie between 0° and Leo “the values 
obtained from the extreme cases where the oxygen does not 


contribute to relaxation and where the oxygen is wholly 


responsible for the relaxation. Due to the short range of the 
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quadrupole interactions, one would expect that o, would lie 
closer to 18° than to 0°. We therefore conclude that (IV 43) 
and (IV 44) will be valid whatever the source of the quadru- 
pole relaxation. 

It is quite generally implied by equations (IV 43) and 
(IV 44) that M1493 must be zero for the ma spin system in 
NaNO,. This is in marked contrast to our experimental results 
which shows that M123 is numerically larger than Mi413° 
Since there is no evidence’ that our crystal is in any way 
defective; we interpret this, surprising.result as ,implying 
that the crystal structure of NaNO, at room temperature is not 
that previously accepted. While our data indicate that the 
By spins are located at sites possessing triad symmetry, the 
orientation of the principal coordinate systems describing the 
quadrupole relaxation is not in agreement with that deduced 
from the crystal structure. Confirmation of. the correctness 
of. the ¢ values is obtained from consideration of the 23Na 
linewidth as a function.of..0. For the. first,setting of.meas- 
urements (at 6 = 2.8°),.the linewidth,was found.to, be. approx- 
imately symmetrical about 9.= 90°. In contrast the second 
set of measurements ($ = 32.8°), the linewidth was found to be 


quite assymetric about 9 = 90°. Comparison with the linewidth 


measurements by Andrew et al. (1962) indicates that ¢ is 
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approximately zero in the first set of measurements and 30° 


in the second set of measurements. 


It*iscconventient)'tomainterpret the M1143 and Mj 193 
values in a slightly different way. We assume that the 


2 ‘ 
ie sites possessetriad symmetry but that the orientation of 


the principal (x) axis is rotated at an angle o* from that 


predicted by the accepted structure. We can therefore write 


2 


re 3 1123 


measured values for Mi4137 —Omdl5etecOs0e8sand M3413 = 04147 + 


roa we find A = 0.187°2-0. 011; and o* ta, 1477372 ier. tne 


cos3¢0 +M sin3¢ in the form Acos3(¢ = $*). Using the 


quantity A may be regarded as a magnitude of the $ -dependent 
term and $* is the rotation necessary to be able to account 
for our data. We note that the angle $~* is remarkably close 
to the value of the angle ¢, considered earlier. Our data 
could at first sight be interpreted by supposing that the Na 
spins are all either type A or type B. However, it is not 
possible to arrange the nitrate groups in. the NaNO. such that 
this can occur. A tentative qualitative explanation of our 
data is that the nitrate groups are rotated from their accepted 
equilibrium positions in such a way that the angle between x 
and x°° is greater than 9, (18°) while the angle between x and 
x°~" is correspondingly. less than to However, measurements 


at the other angles of $ are now required in order to confirm 
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that the ¢ -dependent term really is of the form Acos3(¢ - o*) 
as assumed above. In this connection, measurements at 77 
for » # 0 to complement Niemelas data would be useful. 

Finally we compare the predictions of a very simple 
model with the experimental values of the ratios of the M- 
components. Guided by the accepted configuration of the 
nearest neighbour oxygens, we consider a simple model in which 
the ia nucleus, located at the-orgin of the coordinate 


system, is surrounded by six point charges gq situated at 
(, 7970)» (r,9,2173), > (xr, 0740/3), (©,158,773) 


(571-0) 1)" and” (r,1-0 5773): 


23ya nucleus has therefore triad 


The environment of the 
symmetry and inversion symmetry as is the case in NaNO,. In 
the real crystals the relative motions are very complicated 
and little is known about the relative motion of the eo Na 
nuclei and their nearest neighbours. We therefore make the 
simplifying assumption that the motion of each of the six 
charges relative to the maar nucleus is istropic. In-other 
words, there is no preferential direction in space of the 
relative motions. We now calculate for this simple model the 


fluctuations in the electric field gradient at the origin which 


are caused by this isotropic motion of the point charges. We 
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Caption IV h 


Figure IV h is a schematic representation of the 


charge configuration used in the point charge calculation of 


the various M-components. 
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Shall consider two cases. In the first case we shall calculate 
the changes in the electric field gradient which are linear 
in displacement. In the second case, we shall calculate the 
fluctuation in the electric field gradient which are quadratic 
in displacement. The first case would correspond to the 
anharmonic Raman phonon process and the second case to the 
indirect harmonic Raman phonon process. 

Considering first a simple charge located at (x,y,z), 


we have 


De 3 O) 


V = gq (227-x2-y*) (x7+y742 ) (IV 45) 
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integrated over all relative vibrations. Taking into account 
all six charges and ignoring constants such as e0* in 


equation (II 81) one can write 


Age TO Z 
M3333 = Ey (isan) £" jocv,,) sind do f y dd City 44S) 


where ae Signifies that the expression is summed over the 
configuration of six charges as shown in Figure IV h. After 


some manipulation we find that 


=) Atel eee (1-2c0s70 + 5cos 6). (IV 49) 


We have generalized this expression to the extent of allowing 
Ene vibration an the z-direction to be different from chat in 
hen Xevwand.sdirection. | his.asialblowed, for by introducing a 
factor A for the vibration along the z direction, The value 
X = 1 therefore corresponds to the isotropic case considered 
above. The value A ~*~ 0 corresponds to relative motion pre- 
dominantly in the x-y plane as was suggested by Andrew et al. 
(1962). The value 4 + © represents the situation where the 
relative motion iS predominantly in the z direction along 


the triad axis. Equation (IV 49) becomes 
= 20 cl 2 
M = 18q*u*r "| tse0870-1)? (1-c08%0)-+17 (3-S00e 0) “cos 90 


(IV 50) 


ee hn inh Ry ks Mey a 
a, ‘* YA a i 7 a a 


3eL 


FTAUODOS othe | Pee rT 
nt Sof cs 


(Sk VX) 


eat 1avo. Beronwe at ookesenaie aioe sats soning 
4o3tA' A VE axsip itt ni wots a6, _eogasdo xte, 20.8 
_ tee hag? on noesotuigta 


(ee VI) | 40 "eope + erens-t) * “Super moet <a 
mM aa Dae ps: 


i. 2. 
ae 


oa 


priwolls to tdnasxs eas ast onan Pe 


= 


& pnisuboxsat eed 308 . Bones et outr “snoisoan6. Bas. 
sufsy ent: .agidosyib S aid, al totdsekiv ont x08 a 
bsxeb bedos S260 ra ac 4 | 


=S20 ‘poison ovidatex on 


ent soteiw nositeineas i ala - w 
priols sogso eae s) edt nk 


| O° vos" (o*eone=8)S pate 
toe vr) 


ies 
: loc 


ay i 
7.7 an » & 


a a: a 


4. “a A 


dig) 


Similarly it can be shown for the linear process 


2 


My71= 6a°uer *[sin70[ (3-5sin70) 74 (1/2) (3-5sin20/4) 2+ (3/2) 


(-1+5sin°0/4) ]+i7cos70 [ (5sin“0-1) 7+2(-14+5sin70/4)71] (Iv 51) 


= 0) 
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From symmetry the Mj 42315 zero in this case. We recall that 
the My113 25 measured experimentally is M3 11300839): The 
quantity given in equation (IV 53) is Mo113° AGCOLGdI NG.1Lii-4 ceo 


compare with experiment, this should be multiplied by cos3¢. 


Hence we can write 
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(1-5sin20/4) - (3-5sin20/4) (1-5sin70/4) - (1547/4) (1-5cos70) sin0] 


(IV 54) 
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Figures IV-i, IV 3}, IV kOalidustratei the dependence of Miaii 
M) 3137 M3333 and M1113 on the angle 9 for three values of 

AX = 0.80, 1.00, and 1.25. The large dependence on 0 indicates 
that the local lattice environment has a large effect on the 
M-components. Some indication of the effect of anisotropy 

of. the lattice vibrations is given by comparing the curves 

for different values of i 
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Caption IV 2 


Figure IV i shows the M-components Maui’ M3313, 


M3333 and M1113 as a function of O for the linear process 


and A = 0.80. 
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Caption LV. 4 


Figure IV j shows the M-components Miaau M3137 


and M as a function of 0 for the linear process 
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Caption IV kK 


Figure IV k shows the M-components Miai1° M1313: 


M3333 and M1113 as a function of O for the linear 


process and A = 1.25. 
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2 
X[2(1-5cos“0) + 5sin“0(7cos7@ -1)] + 100cos70sin20(7c0s20-3) 2] 
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Caption Iv: 1 


Figure IV 1 shows the M-components Mijas M1313" 


M3333 and M1113 as a'function of © for the quadratic 


process and A = 0.80. 
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Caption IV m 


Figure IV m shows the M-components Miuil’ Mj 313° 


M3333 and M1413 asa function of © for the’ quadratic case 


and A = 2.00. 
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Caption IV. n 


Figure IV n shows the M-components Miaa1° Mj 3337 


M3333 and M1413 as. a function of-9 for the quadratic case 


and A. = 12125; 
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Since our experiment was concerned with the ratios of 
the M-components and not absolute values of these components, 
and M 


we have plotted M for 


ie Hidde sage imac mie, tian 
he linean case in Figures iV 0, IV p.and IV vq respectively 
tor) =-0.80, 1.00 and .L.25. Similarly,. these, quantities are 
shown, in Figures. IV-r, IVs and IV t for the quadratic case, 
again for A) = 0.80, 1.00 and 1.25. Our experimental values 
are shown on each of graphs IV 0 >.t along with the uncertain- 


ties. In Figure IV q and IV t, we have shown not the ex- 


perimental value of My433/"1111 but the total magnitude of 


the ¢-dependent term defined by Ki? + Me 


ilala le 1123 7 M111" 


is judtifiable since Mj 193 is zero for the theoretical model 


Thass 


we considered. 

We first note that the 0 dependence is similar in the 
linear and quadratic cases for the ratios Mi 3339/1111 and 
M3333 Miq11° However, the ratio M3334 1111 LS uSLOnLPrcantiy, 
smaller in the quadratic case than in the linear case. We see 
that in all cases except for the My373/" 1111 quadratic case, 
the value of 0 which corresponds with the experimentally 
measured ratios is approximately 54° and this happens to be the 
polar angle 90 of the nearest neighbour oxygen atom in the NaNO. 


crystal. The appreciable difference between the theoretical 


linear and quadratic terms for My413/1i11 leads one to wonder 
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Caption IV o 


Figure IV o shows the ratio M as a 


ists” tant 
function of 9 for the linear process.and A = 0.80, 1.00 

and 1.25. Also shown is our experimental value of Mi 313/ 
Miiii ey s oN in NaNO... The shaded region indicates the 


limits of experimental errors. 
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Caption IV p 


Figure IV p shows the ratio M asa 


3333/1121 


function of 90 for the linear process and i = 0.80, 1.00 


and 1.25. Also shown is our experimental value of M3333/ 


Miiil for ena in NaNO... The shaded region indicates the 


limits of the experimental error. 
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Caption IV q 


Figure IV q shows the ratio My33/"1111 asa 
Function of.© forsthe linear process ‘and A = 0.80, 1.00 


and 1.25. Also shown is our experimental value of 


2 Z 23 ‘ 
YM3 413 + MJ 193 y Mii. for the Na in NaNO,. The shaded 


region indicates the limits of experimental error. 
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Caption LV-xr 


Figure IV r shows the ratio M as a 


1313/"1111 


Function of ©: for the quadratic process and 4 = 0.80, 1.00 


and 1.25. Also shown is our experimental value of Mj 313/ 
23 


Mijll Lor Na in NaNO... The shaded region indicates the 


limits of experimental error. 
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Caption IV s 


Figure Iv_s shows the ratio of M asa 


een uae 
function of 0 for the quadratic process and A = 0.80, 1.00 


and 1.25. Also shown is our experimental value of M3 333/ 
23 


Mijii POr Na in NaNO... The shaded region indicates the 


limits of experimental error. 
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Caption.IV t 


Figure IV t shows the ratio Mi 433/"1111 asa function 


of © for the quadratic process and A = 0.80, 1.00.and 1.25. 


, , 2 2 
Also shown is our experimental value of M3113 + Myy93 M1111 


for 2 his in NaNO,. The shaded region indicates the limits 


of the experimental error. 
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whether it is possible to distinguish between a linear and 

a quadratic process on the basis of our crystal measurements. 
Since the agreement between theory and experiment at 54° is 
considerably better in the linear than in the quadratic case, 
it is tempting to conclude that the relaxation mechanism is 
anharmonic. This conclusion must be tempered by the fact 
that the experimental behaviour of the $¢-dependent term is 
quite clear in this case as to the validity of using the 
expression M1137 M3 7130°83¢, with 6 = ican as in developing 
equations (IV 54) and (IV 59). If the 36, term is omitted 
then the agreement is better for the quadratic case than for 
the linear case. 

It is interesting to compare our room temperature data 
results and Niemela's 77°K data in the light of the theoretical 
curves, especially those of Mi 373/ "1111 and M3333/M@)411° Since 
it is known that thermal expansion of NaNO, is some ten times 
greater along the triadaxis than in a direction perpendicular 
to the triad axis (Saini and Mercier, 1934), the polar angle 
of the nearest neighbour oxygen atoms will decrease with in- 
creasing temperature. It follows from Figures IV o, IV p, 

LV x and. IV s that: according to our point seharge: calculation 


the values of My 333/"1111 and M3333/Miq11 should increase with 


the results obtained by Niemela and ourselves. It therefore 
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appears that our point charge calculations are able to provide 
at least a semi-quantitative understanding of our data. 
Finally, relatively small orientation dependence of 
W/W, in NaNo, at room temperature is perhaps not typical of 
the behaviour to be expected in all cases. For example, a 
Ratiowoft M3333/M1701 of near four is expected according to 
our Simple calculation is a crystal where the nearest neighbours 
are located at 9 close to zero. Such a ratio would imply a 
highly anisotropic behaviour of W/W, - An example of a case 
where this may be expected is sodium nitrite (NaNO,). Indeed 
an investigation of the orientation dependence of the relax- 
ation in this crystal would provide a check for the point 
charge model. It would thus appear in retrospect that 
NaNO, was an unfortunate choice of material for the study of 
the orientation dependence of W/W, because the relative dis- 
positions of the nearest neighbour oxygen atoms around the 


sodium sites causes such a small anisotropy in Wo/W)- 


stibvone ed elde sis ancitalaalad oyun dog wo it 
| stab xa0 to piibsstersihay ovises fc oe aso 
(Re Sonsbrsgsh noitetasksto: [lame (lovivelox (yitenkt or "7 
be fesiays don egqsdyeq et owe tersqaigd moon 36 oman: 


a . 

here 
&) yelgqnisxs 1o% sedan ite at betosqxs odes sf _ Te 
oF po Dhuonos bers aie ai thot tear pe sna cael a 


toe 


trodidlp kor jaosat Sit) sxorw Tateyag. 6. at AodtsitoLeo 
\B Yao bivow o1tb= s doue®: .opes a2 geosd =O die swat 
e895 8 to SLlomexs HA - WNW io woiveded, digoxsoaina \y 
heebal -(,0UsW) sd intia me tooe ‘ad bodosgxe od vem aay « 
a | eid to spnebasgep sobtedagt10 oie ‘Se: noisepts 
taieq ot tot dosds ¢ sh bLuow [sdayzo 0 
SoH? Sosqzoxtse) me <804q5 aud bivow 31 ‘<ebom 

Qo rere ont wot Tetsadien. $e sotods ‘etsaudsotn eae 
Veh ovbtsior ond Smcuied r\git ie sonsbnegeb wisn a 
| eat bavers emots nepyxe wwoddp ton sepusbn ods to" ened: 
' a Wyk at Yaomtoe: bas. ima ie: poausa | 


216 


References 


Abragam, A., 1961, The Principles of Nuclear Magnetism 
(Oxford: Clarendon Press), 411. 
Abragam, A., Proctor, W.. G..,11958, Phys.—Rev., 109, 1441. 
Andrew, E. R., 1955, Nuclear Magnetic Resonance (Cambridge 
University Press) 

enarew, Ee RR. ~shades, R. Ge, Hennel ,wJ. Wa, Hughes, D. «Gu, 
L962 puProG. My S SOG..., eV 

AOOGEW, lis, ek, gawansSon ,. k.oM.-, 957 ,.esoc. Phys... Soc.., 70, 436 

BROLEGW ) uli. oR, soworSOn, Ks uM. 1960, Proc. Phys. Soc., 13, D282. 

Bersohn,,R.., 1960,,J3. Chem. Phys., 20, 1505. 

Bloch, .F.., 1946, Bhys. Rev. 10,, 460. 

“Bleembergen, Ny ae urce LL, Ee oi. Pound, Rag. peeoto,) Pays. 
Rev.., g/ 37 19). 

Plume ok. i) aloo ec meneV A OCl as TnSt. 43 7iee) ap 

Edmonds, A. R., 1957, Angular Momentum in Quantum Mechanics 
(Princeton, New Jersey: Princeton University 
Press) 

Goldburg, W. 1., ..959.,. .Bull. sAmen. Phys. SOC 4/5 a ake 

Goldburd, W 1.-, «961, .bnYS..REV = aish 2 el 

Hahn, .E. Le, 949. Phys. Rev., hoe, 145, 


Howling, D. H., 1966, Rev. Sci. Inst., 36, 660. 


| 2 bherioeM mud rou fib muscamot etoRA reer Ses Ka 


“asoneuateat Hs ony sine smith: : 

ad ied sbie ol 

mabsengsM: mene eofgionist ih 410 WA t ae 
ELD |) (ea8e% aobaetreit>/ vrei) vat ak . 

ERRE ROL , . ver agian eel Voor a reas: eae f | 


ae 


Sphbizdin>) sonsnoest oidempeM swetoei .22aer nA ws 

| (sae elesovint nagntys 

Ce ae ee eae eee yeoman: Pe ‘al 

bee , ey , 908 vay’ oot 1SdeL pal bovi pie ite 

GY ..08d /eydt cont (Tees ve = <aoensw2 vet watt 
Jer (1908 cad score ORs a a soatewe \.A a 

26AL OS . age a) tity yOaeL yA 

. Hag aot 14 vot aya ae 

Jeyae ,Sbel (.V uf basen od eecibie 


i =o! = 


ever yee ohaat aioe nhs bs -“ | 
yotexsvinU notesnixd Bice? work smoseonix 8) | 
| (eaert Vy 3 
NEE Vet) CGR Layne som earg <e2er 5 sl ce he Bi 
LE8 SSL, vet -axda oer ¥e% aw 12 
| eb a Von avd seher oh 
083 at deer .to@ «vest E i 6 sent we 

k | 


Za7 


nughes, DD. Gi, 1906," Proc. Phys. Soc. , AS Wired al Hace 

HOgies;. D.9"G., sReea, Ko, L970; Reve scr. Inst. a1 "293 

Hugues) Ds, Gi, Reeds) Ki, Snyder, *R. &., L970, Gane "od - 
Phys., 48, 480. 

Rutsishvili, G. R.;'1956), Publications ofthe Georgian 
TSSeLcure OLesc ences. LV>- 1. 

Niemela, L., 1967, Ann. Acad. Sci. Fennicae A VI, No. 236. 

Ppretala, A., 1963, “Ann..Acad. ‘Sci. Fennicae A. Vi, No.- 2713, 

Pound, R. Toa PE DYS REY. 7/27, 027. 


Pound, Rs 


Vi 

Peund, Re V. 7, 1950, "Phys. Rev. ,; 79, ‘685; 
Vi, Knight, We 2D.) 1950, Rev. Sci. inst. 72172. 
A. 


Redfield, Geo, enys a Rev. 7.98, 1/37. 

Robimeon, Fs) Ni) Ha, 1959, 0. Sci. sInst., 36, 461. 

Saini, ds,emercier, “A.., 1934,) Helv. Phys. Acta. , 7, 267. 

SCnUSteL, (Nien) Loos, Rey...ced  nst., 227 254. 

Snyder, R.-E., Hughes, Ds G., 1970, 0d. Phys. C. (solid 
St. Phys.), in press. 

Van Kranendonk, d.,,1954,) Physica., 20,0782. 

Van Kranendonk, J+; Walker, M.,°1967, Shys. Rev. Letters), 1s, 707. 

Van- Kranendonk;, J.,; Walker, M.; 1968, Can.. J. Phys., 46, 2441. 


Watkins, G. 1D., Pound, Rs Vv.)  1951,. Phys. Rev. 7 62,7) 343. 


ah a ee 
ys 
b> 


Paw oo ay 


E08. £8 \. teat Shoe Lvor over) Sodan Fs pan a8 i: 
U. apd , Over dou A ,sebyae” ax best © 2 a “a 7 
i ue eae ib 

mMipxos) oft to enoissolidsy d80L ee e rh - aye" 

| | Ll ,VE egos Loe 26 ante” | | 
8EL .OUu LEV A epotdnesT .foeg Hae : cn \ Weer vt 


cin Set, BY (A sspianst tod) Been” fet 82eL yh, ye 


.€2@ \8 ,.008 .eyda 2008 ,29eL 1D. ee Ba 


Se 8%) |, em Lay ter iv a 

| 283 ,8e .y we actbadicule ger: ‘ee | 

{O58 ,AS \;taaT . doe! veh ORRE Vy aw Sieh Pee ia ~ 
tl \8¢ oR sais 88@t VO LA st 

£98 \98 , teat . £08 ie PeeE a a ra 

Oa se ligy JA eyag vise | peer ae seloxeu ce x is 

‘BGe Ss . dent -to8 whoa” steet! ay My | 

biter) 2 ayaa Ae ore <“ ia, adppl 

| -eesag ai Cay ie i 

ay OS ,ssdteyae , baer’) ao a ; | 

BL .exosden veal aya \raer om Vaedtete |. t: anobnet 
#i0S \3h oe G m7) 8deE (aM Seniew |, Uae | 
EDE. YS8 v8 soxat Seer ve ame ita: ee 


; ss ree doa - ? 
) 6 The tlre 
“! 


2S 


Whittaker, E., Robinson; G., 1925, The Calculus of 
Observations (Blackie and Sons Ltd.). 

weber, M. d..,. Allen, R. R., 1962, Ampere colloquium, 168. 

Hoekner, EE. <-G., Blumberg, W. Ex, Hahn, %2.-L., 1960, Phys. 
Rev, LIS, 63. 


Mocida, Kk., Moriya, l., £956, 0. Phys. Soc. Japan, 11, 33% 


to agivolsd Prcy seSel (& sot 9 
.(ubdd ea02 bas~eldosts). snotevisedo 
884 murksmpollon sxeqmA Sak ot ey re LLA 
eydd, CdNL , ol .2 Vite » oF ee a 3 é 


Th” Ae, sana Ae faethe 


LED ae om) 


jl ae) s 
i i 
a ; ‘ 
i. , “ ee! 
4 de i ¥ 
i] rf cs ~ 5 % 7 
i y y 
Pres '¢ 
oor ~ 
ween Ve ey 
ee) ; ; 
oh is an a a 4 (ORNS a4 
Dee - 


at hee Lateet 


vet mt Gl" 45). ene ean a 


me's ae view if 


¥ 
> 
ry) 
. 
> 


i ’ om 
- ~ ae 
tint 

ee 


